HLI'HOUUCLU  A I  COVLI'NMLIM  I  l.Xf’LNSL 


KSD-TR-fiS-197 


o 

v— 

00 

00 

to 

T— 

< 

I 

o 

< 


Technical  Report 

724 


Adaptive  Detection  in  Non-Stationary 
Interference,  Part  I  and  Part  II 


t?  ¥ 

xj.j.  ncu^ 


25  June  1985 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lkxim  ton,  Maxsmwsktts 


Prepared  for  the  Deportment  of  the  Air  Force 
under  Electronic  System.  Division  Contract  F 1 9628-85-C-0002. 


Approved  for  Public  rilmr;  distribution  unlimited. 


DTIC 

ELECTE 

AUG  2  91985 


G 


t 


'TrJtPPNSl-- 


The  work  reported  in  this  document  was  performed  at  Lincoln  Laboratory,  a 
center  for  research  operated  by  Massachusetts  Institute  of  Technology,  with  the 
support  of  the  Department  of  the  Air  Force  under  Contract  F1962S  ?5-C-0002. 

This  report  may  be  reproducer!  to  sutisfy  needs  of  U.S.  Government  agencies. 


The  views  and  conclusions  contained  in  this  document  ore  those  of  the  contractor 
and  should  not  he  interpreted  as  necessarily  representing  the  official  policies, 
either  expressed  or  implied,  of  the  United  States  Government. 


The  F.SD  Public  Affairs  Office  has  reviewed  this  report,  and 
it  is  releasable  to  the  National  Technical  Information 
Service,  where  it  will  he  available  to  the  general  public, 
including  foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 

Thomas  J.  Alport,  Major,  USAF 

Chief.  ESD  L  incoln  Laboratory  Project  Office 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


ADAPTIVE  DETECTION  IN  NON  -STATION  ARY 
INTERFERENCE,  PART  I  AND  PART  II 


E.J.  KELLY 
Group  96 


TECHNICAL  REPORT  724 

25  JUNE  19(55 


Approved  for  puhlie  release;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


f\ 

In  this  report,  which  consists  of  two  parts,  the  problem 
of  radar  target  detection  in  a  background  of  non-stationary 
external  interference  is  considered.  The  object  of  the  analysis 
is  to  treat  this  problem  from  the  point  of  view  of  statistical 
decision  theory,  and  to  derive  a  signal  processing  algorithm 
which  accepts  the  totality  of  inputs  on  which  final  decision  is 
to  be  based,  and  performs  both  interference  suppression  and 
target  detection.  It  is  assumed  that  the  radar  is  provided  with 
multiple  RF  input  channels  and  that  target-free  samples,  from 
range  gates  other  than  the  one  in  which  a  target  is  being  sought 
can  be  used  for  the  estimation  of  the  interference  statistics. 

In  Part  I  a  general  formulation  is  given  and  a  likelihood 
ratio  detection  rule  is  derived.  The  probabilities  of  detection 
and  false  alarm  >.re  evaluated  exactly  and  the  performance  of  the 
test  is  illustrated  numerically.  In  Part  II,  a  more  specific 
interference  model  is  introduced,  which  gives  a  more  realistic 
representation  of  the  situations  likely  to  be  encountered  in 
practice,  h  decision  rule  is  derived  which  approximates  the  exac 
likelihood  ratio  tost  for  this  case,  and  approximations  for  the 
defection  <.»nd  false  alarm  probabilities  are  found. 
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PART  I 


(1)  INTRODUCTION 

This  two-part  study  is  devoted  to  certain  aspects  of  the 
problem  of  radar  target  detection  in  a  background  of  external 
interference  which  is  non-stationary  in  character.  In  many 
discussions  of  this  problem,  interference  rejection  is  treated  as 
a  distinct  stage  of  the  processing,  target  detection  taking  place 
as  a  subsequent  operation.  It  is  the  object  of  the  present 
analysis  to  deal  with  this  problem  in  a  unified  way,  leading  to  a 
single  algorithm  which  accepts  the  totality  of  inputs  on  which 
final  decision  is  to  be  based,  and  performs  both  suppression  of 
interference  and  target  detection. 

The  physical  means  with  which  a  radar  copes  with 
interference,  which  usually  take  the  form  of  multiple  RF  input 
channels  and  target-free  data  which  can  be  used  for  the 
estimation  of  the  interference  statistics,  are  taken  as  givens. 
The  target-free  data  is  assumed  to  be  provided  by  the  samples 
from  range  gates  other  than  the  one  in  which  a  taryet  is  being 
sought.  It  is  convenient  to  refer  to  the  former  as  secondary 
data,  while  the  samples  from  the  range  gate  being  processed  for 
targets  are  called  the  primary  data. 

The  analysis  of  Part  1  of  this  study  is  essentially  a 
generalization  of  that  of  the  well-known  paper  ( )  of  Reed, 

Mallett  and  Brennan  (RUB).  These  authors  discuss  an  adaptive 
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procedure  for  the  detection  of  a  signal  of  known  form  in  the 
presence  of  noise  (interference)  which  is  assumed  to  be  Gaussian, 
but  whose  covariance  matrix  is  totally  unknown.  The  possibility 
of  signal  presence  is  accepted  for  the  primary  data,  while  the 
secondary  inputs  are  assumed  to  contain  only  noise,  independent 
of  and  statistically  identical  to  the  noise  components  of  the 
primary  data.  In  the  RMB  procedure,  the  secondary  inputs  are  used 
to  form  an  estimate  of  the  noise  covariance,  from  which  a  weight 
vector  is  determined.  This  weight  vector  i.s  then  applied  to  the 
primary  data  in  the  form  of  a  standard  colored  noise  matched 
filter. 


If  the  output,  of  this  filter  were  compared  to  a  threshold, 
a  complete  detection  procedure  would  bo  obtained.  However,  no 
predetermined  threshold  can  be  assigned  to  achieve  a  given  PFA, 
since  the  detector  is  supposed  to  operate  in  an  interference 
environment  of  unknown  form  and  intensity.  Instead,  in  the  RMB 
paper  an  analysis  is  given  of  the  signal  to  noise  ratio  (SNR)  of 
this  filter  output,  for  given  values  of  the  secondary  data.  This 
SNR  is  a  function  of  the  secondary  data  and  is  therefore  a  random 
variable.  The  probability  density  function  (PDF)  of  this  SNR  is 
deduced,  which  has  the  remarkable  property  of  being  independent 
of  the  actual  noise  covariance  matrix;  it  is  a  function  only  of 
the  dimensional  parameters  of  the  problem. 

In  the  present  study,  the  problem  is  reconsidered  as  an 
exercise  in  hypothesis  testing,  and  the  ad  hoc  RMB  procedure  is 
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replaced  by  a  likelihood  ratio  test.  No  optimality  properties  are 
claimed  for  this  test,  involving  as  it  does  the  maximization  of 
two  likelihood  functions  over  a  set  of  unknown  parameters.  The 
form  of  the  test  is,  however,  reasonable,  and  the  RMB  matched 
filter  output  appears  as  a  portion  of  the  likelihood  ratio 
detection  statistic.  This  test  exhibits  the  desirable  property 
that  its  PFA  is  independent  of  the  covariance  matrix  (level  and 
structure)  of  the  actual  noise  encountered.  This  is  a  general¬ 
ization  of  the  familiar  constant  false  alarm  rate  (CFAk)  behavior 
of  detectors  using  scalar  input  data,  in  which  only  the  level  of 
the  noise  is  unknown.  In  addition,  it  is  shown  that  the  effect  of 
signal  presence  depends  only  on  the  dimensional  parameters  of  the 
problem  and  a  parameter  which  is  the  same  as  the  SNR  of  a 
conventional  colored  noise  matched  filter. 

A  brief  outline  of  Part  I  of  this  study  follows.  The 
detection  problem  is  formulated  in  a  particular  radar  context  in 
Section  2,  where  it  is  also  pointed  out  that  the  actual 
mathematical  analysis  is  considerably  more  general  in  nature.  The 
likelihood  ratio  test  is  derived  in.  Section  3,  where  it  is 
compared  to  the  RMB  procedure  mentioned  above.  In  Section  4,  the 
general  form  of  the  test  is  discussed,  and  its  basic  properties, 
in  particular  its  CFAR  property,  are  exhibited.  The  performance 
of  the  test  is  treated  in  Section  5,  where  the  probabilities  of 
detection  and  false  alarm  are  obtained.  This  performance  is 
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illustrated  numerically  in  Section  6,  which  also  contains  a 
discussion  of  the  results.  Supporting  material,  including  an 
alternate  derivation  of  the  RMB  result,  is  provided  in  four 
appendices. 

The  analysis  in  Part  1  is  self-contained,  dealing  with  a 
single  well-defined  problem.  A  more  specific  (and  more  difficult) 
version  of  the  detection  problem  is  addressed  in  Part  II,  in 
which  the  results  of  Part  I  are  freely  used.  The  analytical 
methods  used  here  lean  heavily  on  the  techniques  of  the  RMB  paper 
(which  in  turn  is  largely  based  on  the  analysis  of  Capon  and 
Goodman^))  with  the  difference  that  the  matrix  transformations 
required  are  carried  out  here  directly  on  the  variables  of  the 
problem,  so  that  much  less  reliance  is  placed  on  known  properties 
of  the  Wishart  distribution. 
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(2) 


FORMULATION  OF  Th If  PROBLEM 


The  mathematical  setting  for  the  formulation  of  this 
detection  problem  will  actually  oe  quite  general,  but  it  is 
introduced  here  first  in  a  relatively  specific  way,  in  order  to 
lend  concreteness  by  way  of  example.  Suppose  that  the  antenna 
system  of  a  radar  provides  a  number,  say  M,  of  RF  signals.  These 
may  be  the  outputs  of  array  elements,  subarrays,  beamformers  or 
any  mix  of  the  above.  The  radar  waveform  is  supposed  to  be  a 
simple  burst  of  identical  pulses,  say  j  in  number,  and  target 
detection  is  to  be  based  upon  the  returns  from  this  burst. 

Further  processing  is  possible,  of  course,  and  the  'detections' 
of  the  burst  processor  may  be  taken  as  inputs  for  subsequent 
binary  integration  over  a  string  of  bursts.  In  any  case,  the 
burst  processor  makes  a  decision,  comparing  some  function  of  its 
input  data  (called  the  detection  statistic)  to  a  threshold,  and 
the  design  of  this  processor  is  the  same  as  if  its  decision  were 
to  be  the  final  one. 

In  effect,  the  radar  front  end  carries  out  amplification, 
filtering  and  reduction  to  base  band,  at  which  point  the 
quadrature  signals  are  subjected  to  pulse  compression,  the  final 
stage  of  filtering.  The  order  in  which  these  things  are  carried 
out  is  immaterial  to  the  present  model,  since  rt  is  not  addressed 
to  the  problems  of  realization  and  channel  matching,  although 
these  are  of  great  importance  in  practice.  The  in-phase  and 
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quadrature  output  pairs  are  next  sampled  to  form  range  gate 
samples  for  each  pulse,  say,  G  range  gates  per  pulse  .  This 
results  in  a  total  of  MJG  complex  samples  for  the  burst.  Signal 
presence  is  sought  in  one  range  gate  at  a  time,  hence  the  primary 
data  consists  of  the  MJ  samples  from  a  single,  unnamed  range 
gate.  These  samples  are  arranged  in  a  column  vector,  z,  of 
dimension  N  =  MJ  ■  The  secondary  data  consist  of  the  outputs  of  K 
range  gates,  forming  a  subset  of  the  G-l  remaining  ones,  and 
these  are  described  by  the  set  of  vectors,  z(k),  (k  =  1...K).  The 
decision  rule  will  be  formulated  in  terms  of  the  totality  of 
input  data,  without  the  a  priori  assignment  of  different 
functions  to  the  primary  and  secondary  inputs. 

The  secondary  data  are  assumed  to  be  free  of  signal 
components,  at  least  in  the  design  of  the  algorithm,  and  any 
selection  rules  applied  to  make  this  assumption  more  plausible 
are  ignored.  The  primary  data  may  contain  .«  signal  vector, 
written  in  the  form  bs ,  where  b  is  an  unknown  complex  scalar 
amplitude,  and  s  is  a  column  vector  of  N  components  describing 
the  signal  which  is  sought.  The  modeled  variation  of  signal 
amplitude  and  phase  among  the  array  inputs  is  included  in  s,  as 
well  as  pulse  to  pulse  variations,  such  as  those  relating  to  a 
particular  target  doppler  velocity.  The  problem  of  unknown 
doppler,  or  other  unknown  signal  parameters,  is  mentioned  briefly 
below.  It  should  be  noted  that  the  signal  vector,  s,  can  be 


normalized  in  any  convenient  way,  since  an  unknown  amplitude 
factor  is  already  included,  and  we  retain  the  freedom  to  assign  a 
norm  to  s  at  a  later  point,  where  it  will  be  most  advantageous  to 
make  a  specific  choice. 

The  total  noise  components  of  the  data  vectors, 
representing  all  sources  of  internal  and  external  noise  and 
interference,  are  modeled  as  zero-mean  complex  Gaussian  random 
vectors.  The  noise  component  of  the  primary  vector,  z,  is 
chc  acterized  by  the  unknown  covariance  matrix,  M.  Each  of  the 
z(k)  is  assumed  to  share  this  NxN  covariance  matrix,  and  the 
vectors  z  and  the  z(k)  are  all  mutually  independent.  All  Gaussian 
vectors  are  assumed  to  have  the  'circular'  property  usually 
associated  with  I  and  Q  pairs. 

The  key  features  of  this  model  are  the  Gaussian 
assumption,  the  independence  of  the  primary  and  secondary  inputs, 
and  the  assumption  that  these  share  a  common  covariance  matrix. 
The  structure  of  the  N-vectors,  in  particular  the  doubly  indexed 
model  used  to  describe  multiple  pulses  and  multiple  array 
outputs,  is  not  used  in  the  following.  One  may  equally  well  think 
of  J  as  being  to  unity,  in  which  case  the  analysis  corresponds  to 
a  situation  in  which  detection  is  based  on  the  returns  from  a 
single  pulse. 
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(3) 


THE  LIKELIHOOD  RaTIO  TEST 


Consider  a  single  input  vector  from  the  secondary  data 
set,  say  z(k).  If  the  covariance  matrix  of  this  vector  is  M: 

M  =  Ej  z(k)  z(k)f } , 

then  the  N-dimens ional  Gaussian  PDF  of  this  complex  random  vector 
will  be 

,Wk)]=  i_  •-  *«' . 

In  the  notation  used  here  the  double  bars  signify  the  determinant 
of  a  matrix,  and  the  superscript  dagger  symbolizes  the  conjugate 
transpose  of  a  vector.  Each  of  the  secondary  data  vectors  has 
this  same  PDF,  and  under  the  'noise-alone'  hypothesis,  the 
primary  vector  does  so  as  well,  hence  the  joint  PDF  of  all  the 
input  data  is  the  product: 

fo[z,2(l)»",2(K)]  =  f[2]  1 1  f[z0<)]  • 

If  v  is  any  N-vector,  we  can  write  the  following  inner 
product  in  the  form  of  a  matrix  trace  (Tr): 

vfMHv  =  Tr(M^V) 

where  V  is  the  open  product  matrix 

V  =  vvt 

When  this  equivalence  is  applied  to  all  the  factors  of  the 
joint  PDF,  it  will  be  seen  that  the  latter  may  be  written  in  the 
convenient  form 
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where 


f0l>*zO) 


1  «-  Tr(M_1  T0) 

ttn  1!  m  I! 


K+1 


=_JL/i2UV  :(t)i(k)’] 

K+1  \  fei  > 


Under  the  1  signal-plus-noise '  hypothesis,  the  z(k)  have 
the  same  PDF  as  before,  and  the  PDF  of  the  primary  vector  is 
obtained  by  replacing  z  by 


-  z-bs. 

The  resulting  joint  PDF  of  the  inputs  is  then 


where  now 


f1[z,z(l)...,z(K)] 

Ti  =  kVi  ( (7-b3)^bs)t  + 


e-  Tr(M_1  T^j  K+1 


In  the  likelihood  ratio  testing  procedure,  the  PDF  of  the 
inputs  is  maximized  over  all  unknown  parameters,  separately  for 
each  of  the  two  hypotheses.  The  ratio  of  these  maxima  is 
the  detection  statistic,  and  the  hypothesis  whose  PDF  is  in  the 
numerator  is  accepted  as  true  if  it  exceeds  some  preassigned 
threshold.  The  maximizing  parameter  values  are,  by  definition, 
the  maximum  likelihood  (MI.)  estimators  of  these  parameters,  hence 
the  maximized  PDF's  are  obtained  by  replacing  the  unknown 


parameters  by  their  MI.  estimators. 
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We  begin  with  the  noise-alone  hypothesis,  maximizing  over 
the  unknown  covariance  matrix,  M.  Of  all  positive  definite  M 
matrices,  the  one  which  maximizes  the  expression  inside  the  curly 
brackets  of  this  PDF  is  simply  Tq.  This  in  equivalent  to  the 
statement  that  the  ML  estimator  of  a  covariance  matrix  is  equal 
to  the  sample  covariance  matrix,  which  is  well  known(3>.  When 
this  estimator  is  substituted  in  the  PDF,  the  trace  which  appears 
there  becomes  the  trace  of  the  NxN  unit  matrix,  which  is  just  N, 
and  we  find 

/  j  \  K+1 

MmX  *°  \  (®tr)M  ||  T0  |  /  ' 

The  same  procedure,  applied  to  the  siynal-plus-noisc 
hypothesis  yields  the  formula 

/  1  \ K+1 

Msx  ■ 

and  it  remains  to  maximize  this  expression  over  the  complex 
unknown  signal  amplitude,  b.  Since  b  appears  only  in  this  PDF,  we 
can  form  a  likelihood  ratio,  L(b),  at  this  point  and  subsequently 
maximize*  it  over  b.  It  is  more  convenient  to  work  with  the 

(K+l)st  root  of  this  ratio,  arid  we  put 

K+1 

L(b)  «  U(b)| 

Obv icusly , 


=  LLlII 
II  T, 
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and  the  final  likelihood  ratio  test  takes  the  form 


Max  4(b)  = 
b 


II  Tn  || 

Min  ||  Tj  || 
0 


>  l 


o  • 


The  threshold  parameter  on  the  right  will  evidently  be  greater 
fhan  unity,  since  the  denominator  on  the  left  equals  the 
numerator  for  the  choice  b  =  0,  and  we  are  maximizing  over  b. 
To  proceed,  we  define  the  matrix 

K 

S  =  £  i(k)  x(k)» 
k=1 


which  involves  only  the  secondary  data.  This  macrix  is  K  times 
the  sample  covariance  matrix  of  these  data,  and  it  satisfies  the 
well-known  Wishart  distribution.  The  only  property  of  this 
distribution  that  we  need  here  is  the  fact  that  for  K  >  N,  a 
condition  we  now  impose,  the  matrix  s  is  non-singular  with 
probability  one.  S  is,  o.t  course,  positive  definite,  and  hence 
Hormitian.  Wo  use  a  lemma  proved  in  Appendix  B  to  evaluate  the 
determinants  of  both  sides  of  the  equation 

(K+1)  T0  =  S  +  zzt 

with  the  result 

(K+1)N  II  To  II  =  II  s  II  (l  +  zV'z). 

Similarly,  wo  have 

(K+1)N  ||  Tt  ||  =  ||  S  ||  (l  t  (z-bs)tS-1(z-bs)j  . 

Now  is  a  good  time  to  minimize  this  quantity  over  b,  and  wc 
do  this  by  comp le Liny  the  square: 
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is  equivalent  to  the  test 


We  note  that  tiq  lies  between  the  values  zero  and  one. 

If  the  target  model  is  generalized,  so  that  the  signal 
vector  still  contains  one  or  more  unknown  parameters  (such  as 


target  dopplor),  the  likelihood  ratio  obtained  above  must  next  be 
maximized  over  these  parameters.  It  is  clear  that  this  is 
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equivalent  to  maximizing  ri  itself  over  the  remaining  target 
parameters.  This  maximization  generally  cannot  be  carried  out 
explicitly,  and  the  standard  technique  is  to  approximate  it  by 
evaluating  the  test  statistic,  in  this  case  n ,  for  a  discrete  set 
of  target  parameters,  forming  a  'filter  bank',  and  declaring 
target  presence  if  any  filter  output  exceeds  the  threshold.  Our 
purpose  in  discussing  this  here  is  only  to  show  how  our  test  can 
be  generalized  in  this  straightforward  way,  but  from  now  on  we 
ignore  any  additional  target  parameters,  which  is  equivalent  to 
concentrating  on  the  performance  of  a  single  member  of  the  filter 
bank. 

A 

For  comparison  with  the  RMB  procedure,  we  introduce  M,  the 
ML  estimator  of  the  noise  covariance,  based  on  the  secondary  data 
alone.  We  have  already  noted  that  this  estimator  is  equal  to 

5=ls. 

The  likelihood  ratio  test  can  then  be  written  in  the  form 

|(stM"tz)|2 

—r—. - ; — 7^-7—  >  • 

(s*M-Is)  [1  + 

We  note  that  the  secondary  inputs  enter  this  test  only  through 

A 

the  sample  covariance  matrix,  M,  and  also  that 

(s^M"4*)  -  (wt  1) 

where  w  is  the  RMB  weight  vector 
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W  a  M  1  S  . 


The  RMB  test  itself  is  just 

|(w  z)| 8  >  threshold  „ 

which  has  the  form  of  the  colored  noise  matched  filter  test,  with 

A 

M  replacing  the  usual  known  covariance  matrix  of  the  noise. 

The  presence  of  the  signal-dependent  factor  in  the 
denominator  of  the  expression  for  n  causes  this  detection 
statistic  to  be  unchanged  if  the  signal  vector  is  altered  by  a 
scalar  factor.  Since  the  normalization  of  this  vector  has  been 
left  arbitrary,  this  invariance  is  highly  desirable.  In  effect, 
this  factor  in  the  denominator  is  normalizing  s  for  us,  in  terms 
of  the  estimated  noise  covariance.  The  entire  detection  statistic 
is  also  invariant  to  a  common  change  of  scale  of  all  the  input 
data  vectors,  a  minimal  CFAR  requirement.  Further  properties  of  n 
will  be  developed  in  the  following  section. 

In  the  limit  of  very  large  K,  one  expects  the  estimator, 

A 

M,  to  converge  to  the  true  covariance  matrix,  M,  at  least  in 
probability.  Moreover,  it  can  be  shown  that  the  quantity 

(z'M-'l) 

an  inner  product  utilizing  the  actual  covariance  matrix  instead 
of  its  estimator,  obeys  the  chi-squared  distribution,  with  2N 
degrees  of  freedom,  and  hence  this  term,  when  divided  by  K, 
converges  to  zero  in  probability,  as  K  grows  without  bound.  In 
this  sense  the  likelihood  ratio  test  passes  over  into  the 
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conventional  colored  noise  matched  filter  test,  as  the  number  of 
sample  vectors  in  the  secondary  data  set  becomes  very  large. 
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(4) 


PROPERTIES  OF  THE  LIKELIHOOD  RATIO  TEST 


The  likelihood  ratio  test  will  be  discussed  in  terms  of 
th  vandom  variable  n,  the  decision  statistic  eventually  obtained 
in  >. ?  preceding  section.  The  definition  of  n,  as  well  as  that  of 
tht  matrix  S  on  which  it  depends ,  are  reproduced  here  for 
coi  \  •'  1  i  ence; 

l(»,s-‘z)l2 

”  (.V1.)  [1  +  (z's-’z)] 

K 

s  =  Yl  z0<)  z(k)*  • 

k=1 

The  random  variable  n  is,  of  course,  a  function  of  both  the 
primary  and  secondary  data,  and  as  a  preliminary  to  discussing 
its  actual  PDF,  some  useful  properties  are  first  derived.  We 
begin  with  the  noise-alone  case,  and  assume  that  the  actual  noise 
covariance  matrix  is  M. 

The  matrix  M  is  positive  definite,  and  hence  a  positive 
definite  square  root  matrix  can  be  defined.  Since  M  can  be 
diagonalized  by  a  unitary  transformation,  it  can  be  represented 
in  the  form 

M  =  U  A  U* 

where  the  columns  of  the  unitary  matrix,  U,  are  eigenvectors  of 
M,  and  A  is  diagonal.  The  diagonal  elements  of  A,  say  Mn), 


(n=l...N),  are  the  real,  positive  eigenvalues  of  M.  In  case  of 
degeneracy  of  an  eigenvalue,  the  corresponding  eigenvectors  are 
assumed  to  have  been  orthogonalized .  The  square  root  may  be 
defined  by  the  representation 

=  uaM 

where  A^/2  is  diagonal,  with  diagonal  elements  [l(n)  ]^/2. 

M-1/2  is  similarly  defined  in  terms  of  A“l/2,  and  it  is 
easily  seen  to  be  the  inverse  of  M^-/2.  Uniqueness  of  the  square 
roots  is  not  necessary  for  our  purpose,  only  their  existence  and 
positive  definite  (hence  also  Hermitian)  character. 

Now  consider  the  vector 


and  the  similarly  transformed  secondaries 

f(k)  a  M-^  z(k) . 

The  new  vectors  are  zero-mean  Gaussian  variables,  but  with 
covariance  matrix  equal  to  1^ ,  the  NxN  identity  matrix.  This 
follows  directly  from  the  definitions: 

E$zzf|  M“&  =  M"fc  M  =  IN  . 

with  identical  reasoning  for  the  transformed  secondaries.  The 
linear  transformation  introduced  here  is,  of  course,  a  whitening 
transformation. 

We  note  that  the  scalar,  n,  depends  on  the  data  and  signal 
vectors  only  through  inner  products.  By  inverting  the  whitening 
transformation  we  may  evaluate,  for  example,  the  product 


(ztS“1z)  =  (ff  r)  =  (V  (M'isM'i)-1  . 

We  define  the  new  matrix 

J  m  S 

and  substitute  for  S,  finding 

K  t  .  K 

y  =  Y'  M~S  z(k)  z(k)*  M~2  =  Yj  f(k)  ^(k)f  . 
k=1  k=1 

Therefore,  the  new  S-matrix  is  K  times  the  sample  covariance 
matrix  of  the  whitened  secondaries,  and  the  random  variable 

E  ^  (ztS~lz)=  (?V~V) 

is  seen  to  be  independent  of  M,  being  expressible  as  a  function 
of  K+l  independent  Gaussian  vectors,  each  of  dimension  N,  and 
each  sharing  the  covariance  matrix,  In*  The  PDF  of  l,  like  that 
of  the  RMB  signal-to-noise  ratio,  is  therefore  a  universal 
function  of  the  dimensional  parameters,  N  and  K,  alone. 

The  other  inner  products  in  the  decision  statistic 
are  handled  in  an  analogous  manner;  thus 

(•fs_lz)  =  (sViyV)  =  (<V~V) 

where  t  stands  for  the  whitened  signal  vector 

At  this  point  we  make  the  deferred  definition  of  signal 
normalization,  by  taking  t  to  be  a  unit  vector: 

(<f/)  =  (sVs)  H  1 . 

This  choice  gives  specific  meaning  to  the  signal  amplitude 
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parameter,  b,  whose  square  is  now  a  proper  signal-to-noise  ratio, 
for  which  we  introduce  the  symbol  a: 


-|b|*  =  (EWtM-‘Ew). 


When  the  obvious  substitutions  are  made  in  the  final  inner 
product,  we  obtain 

”  «V-‘ 0  [1  +  ' 

The  dependence  on  M  is  now  confined  to  t,  and  it  will  be  shown 
below  that  even  this  dependence  on  the  true  covariance  matrix  is 
illusory.  When  a  signal  is  present,  z  is  replaced  by 

z  =  bs  +  n 

where  n  has  all  the  properties  attributed  to  z  in  the  noise-alone 
case.  In  this  situation,  the  whitened  data  vector  is 


whe  re 


y  =  M  kz^bt  +  v 
v  =  n  , 


which  is  statistically  identical  to  the  whitened  data  vector  in 
the  noise-alone  case.  We  have  therefore  found  that  when  signal  is 
present,  the  PDF  of  n  depends  on  M  only  through  b  and  t,  and  the 
dependence  on  the  unit  vector,  t,  is  again  only  apparent,  as  we 
now  proceed  to  show. 

Suppose  the  whitening  transformation  is  followed  by  a 
unitary  one,  in  which  the  whitened  vectors  are  expressed  as  the 
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products  of  a  unitary  matrix  and  a  new  set  of  random  vectors. 
These  new  random  vectors  are  statistically  indistinguishable  from 
their  predecessors,  and  it  would  only  be  confusing  to  introduce  a 
new  notation  for  them.  Tracing  this  transformation  through  the 
inner  products,  we  find  that  only  the  normalized  signal  vector  is 
changed:  t  is  replaced  by 

b  ly 

where  is  the  unitary  matrix  characterizing  this  last 
transformation.  Any  unit  vector  in  the  complex  M-space  can  be 
realized,  as  ti ,  by  such  a  transformation.  In  particular,  we 
can  cause  ti  to  be  a  'coordinate  vector',  for  which  a  single 
element  is  unity,  the  remaining  (N-l)  elements  vanishing.  It  is 
for  this  reason  that  the  PDF  of  n  depends  on  M  only  through  the 
meaning  of  the  signal  amplitude  parameter,  b.  In  fact,  this  PDF 
can  depend  only  on  b,  N  and  K,  and  hence  the  false  alarm 
probability  of  the  likelihood  ratio  detector,  namely 

PFA  =  Prob$77  >  7j0| , 

is  independent  of  M,  and  this  is  the  generalized  CFAR  property 
claimed  in  section  1, 
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(5)  PROBABILITY  DISTRIBUTION  OF  THE  TEST  STATISTIC 

We  now  take  advantage  of  our  freedom  to  make  a  unitary 
transformation,  and  choose  for  tj.  a  vector  whose  first  element 
is  unity,  all  others  being  zero.  This  can  be  accomplished  by 
choosing  for  U]_  a  matrix  whose  first  row  is  the  conjugate 
transpose  of  t,  and  whose  other  rows  are  the  conjugates  of  unit 
vectors  orthogonal  to  t.  Understanding  that  this  choice  has  been 
made,  we  drop  the  subscript  on  t,  so  that  n  is  still  given  by  the 
formula  of  the  preceding  section. 

This  form  for  t  makes  it  expedient  to  decompose  all 
vectors  into  two  components,  an  A-component  consisting  of  the 
first  element  only,  and  a  B-component  consisting  of  the  rest  of 
the  vector.  Thus  we  write 

-IS] 

where  the  A-component  is  a  scalar  and  the  B-component  is  an 
(N--1 ) -vector .  In  this  notation,  the  signal  vector  is  just 


the  zero  being  (N-l )  dimensional.  Matrices  are  decomposed  in 
analogous  fashion,  and  we  write 


Note  that  the  AA-element  is  a  scalar,  the  BA-element  is  an  (N-l) 
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dimensional  column  vector,  and  so  on.  We  also  give  a  name  to  th 
inverse  of  this  matrix,  decomposing  it  as  well: 


-w  =  [ 


^AA  ^AB 
^BA  ^BB 


With  this  notation  we  have,  simply, 

(Stf-U)  =  (S?t)  =  . 


while 


«V-V)  =  [io]f 

-  •'ba  •'bb  J  l  fa 


It  is  important  to  keep  in  mind  that  we  now  have  a  four-fold 
decomposition  of  the  total  input  data  set  into  primary  and 
secondary  vectors,  each  of  which  is  divided  into  A-  and 
B-components . 

According  to  the  Frobenius  relations  for  partitioned 
matrices , 

■^AA  =  (*^AA  ”  "Aa )  * 

which  is  a  scalar,  and  also 

*.a  =  • 

Since  the  NxN  sample  covariance  matrix  and  its  inverse  are 
Hermitian,  we  obtain 
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and  therefore 


At 


_  t  =  _ 
•'•a 


^aa^ab^w 


-l 


V)  =  ^AA^A  "  Wm  V.)  . 


The  ftnal  inner  product  is  expanded  as  follows: 

(/y-1,)  =  *A>/  +  2  . 

where  we  have  applied  the  identity 


(ufv)  =  (v^u)* 

to  the  (N-l )  dimensional  inner  product 

r.'rK 

Next,  we  complete  the  square  in  this  last  expression,  writing 


(^y  f)  -  1 1 k  +  yAA  Ja»?»I 

+  —  ^AA  ^BA^Ab)?- 1  • 

and  by  using  a  Frobenius  relation  in  reverse  we  see  that 

<P  _  €>  e>  <P  — 

*^AA  J»A*rA»  ~  *  •§ 

Finally,  combining  these  results,  we  find 

=  ,M 1,4  -  -WV  +  (r,V»,)  • 

When  those  evaluations  are  substituted  into  our  expression 

for  n,  the  result  can  be  expressed  in  the  apparently  simple  form: 

„_X__ 

v  1  +  X  +  S, 

Wo  have  introduced  hero  the  notation 

=  Vs)  1 

which  will  be  retained,  and  the  temporary  notation 
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^  ~  ?»l  * 

Note  that  %  is  just  like  the  quantity  E  defined  earlier, 
except  that  the  dimensionality  of  the  vectors  involved  is  now 
N-l .  The  last  form  of  the  decision  test,  namely 

V  >  • 

is  evidently  equivalent  to 

x>-j^Vo<1 

We  shall  leave  the  test  in  this  form  for  a  time,  while  we  examine 
the  statistical  properties  of  the  quantities  which  enter  into 
it. 

A  previous  evaluation  for  the  leading  factor  in  X  can  be 
used  to  obtain  the  following  form: 

„  \tk  ~  ^AB^BB  ^>1 

x  =  - — rr -1  :r 

aAA  “  'yAB'/6B  ^BA 

Wo  make  use  of  the  definitions  to  express  the  denominator  as  a 
sum: 

^AA  “  '^AH^BB  "^BA  =  ^  ^/aM  ”  ’-'ab'^BB  ^  aM  * 

This  is  the  same  as  the  sum  of  squares 

K 

because  the  terms  supplied  to  complete  this  square  add  up  to 


ze  ro : 


jjp  ^A(k)  ^AB'^BB  ?s(k)j 

K  K 

=  E  *A(k)»’B(*Ot‘'w”1^tt  -  ^ab^bb'1  E  ?B(kWk)f  ^b'^ba 

k=1  k=1 


•^AB^BB  ^BA 


■^AB^BB  "^BB^BB  ^BA 


=  0  . 


The  evaluation  of  the  suras  here  follows  from  the  definitions  of 
the  partitioned  matrix  elements.  We  introduce  the  notation 

y(k)  =  ?A(k)  -  ■^Ai^ii  fi(k) 

for  the  terms  of  the  sum,  and  the  analogous  notation 

y  b  ^a  _  ^a*^»» 

for  the  quantity  appearing  in  the  numerator  of  X,  so  that  the 
likelihood  ratio  test  can  be  written  in  the  more  explicit  form 


y-  _!£_ 

I  l/(k)!2 

k-1 


> 

1  _  Vo 


(l  +  £.)  • 


We  proceed  by  fixing  the  b- vectors  temporarily,  and 
consider  the  probability  densities  of  all  quantities  entering 
into  the  decision  statistic  to  bo  conditioned  on  these  values. 
The  conditional  probabilities  of  detection  and  false  alarm  will 
be  evaluated  first,  and  the  condition  will  then  be  removed  by 
taking  expectation  values  over  the  joint  1'bF  of  the  b-vectors. 
With  the  b-vectors  fixed,  only  the  K-t  1  scalar  A-compononts  are 
random,  and  we*  show  now,  under  this  condition,  that  y  and  the 
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y(k)  are  Gaussian  variables,  that  y  is  uncorrelated  with  the 
y ( k ) ,  and  that  the  latter  have  a  covariance  matrix  with  simple 
properties . 

Using  the  definitions  of  the  y's  and  of  the  AB  matrix 
element  which  enters  there,  we  can  express  these  quantities  in 
the  form 

K 

y =  ~  "Yj  T'aMt'bOO^b* 

k=1 

and 

K 

y(k)  =  *A(k)  -  g  nO)^(i)VIB“V((k)  . 

This  represents  the  y's  as  linear  comb ina lions  of  the 
A-componcnts ,  and  hence  proves  their  conditional  Gaussian 
character.  Moreover,  the  y(k)  have  zero  mean  in  all  cases,  while 
the  conditional  mean,  written  bn,  of  y  in  the  general  case  is 

y  —  ^  t  k  —  b  • 

as  a  result  of  our  choice  of  signal  vector. 

The  linear  dependence  of  the  y's  on  the  A -components  is 
best  expressed  in  terms  of  the  quantities 

q(k)  a 

and 

Q(i.k)  *  rB( OV.fV.W  • 

which  are  constants  under  the  conditioning.  Obviously, 

K 

y  =  **  -  E  **00 

k=1 


1-26 


and 


K 

y(t<)  =  fA0<)  - 1]  **0)  Q0>k)  • 

k=1 

The  q(k)  may  be  considered  as  the  components  of  a  K-vector,  q, 
and  the  Q(i,k)  as  the  elements  of  a  KxK  matrix,  Q.  The  desired 
properties  of  the  y's  flow  from  the  following  facts  about  this 
new  vector  and  matrix: 

Qq  =  q 

and 

Q2  =  Q  . 

To  prove  the  first  of  these,  we  write  it  out  in  component 

form: 


K  K 

JC  Q(k,i)q(0  •-=  g  iWO*  V1*.  • 


The  sum  over  i  regenerates  the  BB  matrix  element: 

K 


*(»)*(»)*  =  . 


[Hi 


(as  happened  when  the  denominator  of  X  was  expresed  as  a  sum  of 
squares),  and  the  result  follows  immediately.  The  idempotent 
character  of  Q  is  proved  in  the  same  way.  We  also  note  that  Q  is 
Hormitian,  and  that  its  trace  is  N-l: 


K 


Tr(Q)«  g^,(k)T^M">t(k) 

/  K 

=  Tr (v1  5]  *,(10*00* 

'  k-1 

=  Tr(IN  !>  =  N  "  1  • 
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Note  that  we  are  dealing  with  the  trace  of  a  KxK  matrix  on  the 
left  side  here,  and  of  (N-l)x(N-l)  matrices  on  the  right.  All  of 
these  results  will  be  required  in  the  following. 


and  hence 

E,|y-bf  =  l  +  S.  . 

This  last  result  is  responsible  for  a  significant  simplification 
of  the  statistics  of  the  likelihood  ratio  test. 


1-28 


Finally,  we  compute  the  conditional  covariance  of  the 
y(k).  We  use  the  notation  <5(i,k)  for  the  elements  of  the  unit 
matrix,  so  that  y(k)  can  be  written 

K 

y(k)  =  E  r*(0  [«('*)  "QC*)]  . 

Using  the  independence  of  the  A-variables  again,  we  obtain 

K 

Ey(k)y(n)*  =  [<5(l,k)  -  Q(l,k)]  [«5(ltn)  -  Q(l.n)]* 

J-d 

K 

=  <5(k,n)  -  Q(n,k)  -  Q(k,n)*  +  V  Q(!.k)Q(l,n)*  . 

1=1 

Since  Q  is  Hermitian  and  idempotent,  we  find  the  simple  result 

Ey(k)y(rt)*  =  <5(n,k)  -  Q(n,k)  . 

The  likelihood  ratio  test  is  now  rearranged  slightly  to 

read 


-^L  >_5°-  y  |y(k)|2. 

In  view  of  fact  that  the  conditional  variance  of  y  equals  the 
denominator  on  the  left,  it  makes  sense  to  define  a  normalized 
variable 

Conditioned  on  the  B-vectors,  w  is  Gaussian  and  independent  of 
the  y(k).  It  has  a  conditional  variance  of  unity,  and  a 
conditional  mean  value: 
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ElW  =  (i  +  s1)i7i  ' 

In  the  noise-alone  case,  the  conditioning  has  no  effect  on  the 
PDF  of  w.  The  sum  over  k  is  also  given  a  name: 

K 

t  -  E  i/oof  • 

k=1 

and  the  test  is  now  written 

|wjZ>(/0-l)T  * 

where  the  original  threshold  constant  has  been  reintroduced.  In 
fact,  it  is  easily  verified  that  our  original  likelihood  ratio 
given  by 

w2 

l  =  1  +—  . 

We  now  turn  to  the  properties  of  T.  Given  the  B-vectors, 
the  joint  PDF  of  the  y(k)  is  zero-mean  Gaussian  with  eova;.  ince 
matrix  J: 

J(l.k)  a  <5(1, k)  -  Q(k,l)  . 

The  conditional  characteristic  function  of  T  is  therefore 

«,(X)  ^  E.Je,XT|  =  ||I  -  lXJf1  . 

Since  Q  is  idompotent,  its  eigenvalues  are  either  zero  or 
one,  and  from  the  value  of  its  trace  we  see  that  Q  must  have 
exactly  N-l  unit  eigenvectors.  It  follows  that  >1  has  K+l-N  unit 
eigenvectors,  the  others  being  zero,  and  thus 

4,(X)  =  (i  -  IX)-(K+1“N> 

This  is  the  characteristic  function  of  a  chi-squared  random 


is 
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variable/  and  the  PDE'  of  T  is  simply 


rK-N 


*.(T)  = 


— T 


(K-N)l 


It  is  remarkable  that  the  statistical  properties  of  T  are 
independent  of  the  actual  values  of  the  conditioning  B-vectors, 
and  we  can  consequently  drop  the  subscript  on  its  PDF.  Moreover, 
T  is  statistically  equivalent  to  the  sum  of  the  squares  of  K+l-N 
independent,  complex  Gaussian  variables,  each  of  which  has  zero 
mean  and  unit  variance.  If  we  let  w(k),  ( k=l ... K+l-N ) ,  be  such  a 

set,  then  T  is  statistically  indistinguishable  from  the  sum 


K+l-N 


|w(k)|a  . 


The  properties  of  the  likelihood  ratio  test  are  therefore 
identical  to  the  properties  of  the  simple  test 

K+l-N 

IwfX'o-DE  1*001*  • 

where  the  w(k)  are  now  also  taken  to  be  independent  of  w.  Tne 
probability  of  the  truth  of  this  inequality  is  still  conditioned 
on  the  B-vectors,  but  this  conditioning  appears  only  through  the 
quantity  Ejj,  which  is  contained  in  the  conditional  mean  of  w. 

This  equivalent  decision  rule  represents  the  behavior  of  a 
simple  scalar  CFAR  test,  in  which  the  power  in  one  complex  sample 
(a  single  radar  hit),  being  tested  for  signal  presence,  is 
compared  to  a  threshold  proportional  to  the  sum  of  the  powers  of 
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K+l-N  samples  of  noise.  This  problem  is  quite  familiar,  and  the 
test  just  described  is  also  a  likelihood  ratio  test  in  the 
corresponding  situation.  The  performance  of  the  scalar  CFAR 
detector  is  very  simple,  and  in  particular,  its  PFA  is  just 


(ir  - 


In  this  case,  when  the  signal  amplitude  is  zero,  the  conditioning 
B-vectors  do  not  appear  at  all,  and  hence  this  simple  formula 
gives  the  PFA  for  our  original  likelihood  ratio  test. 

The  probability  of  detection  (PD)  of  the  scalar  CFAR  test 
is  also  well  known,  and  in  our  case  it  depends  on  the  conditional 
SNR,  which  is  the  squared  magnitude  of  the  conditional  mean  of 
w.  In  terms  of  the  colored  noise  matched  filter  SNR,  a,  defined 
earlier,  and  the  quantity 

r  =  — * 


1  +  E. 


this  conditional  SNR  is  just  ra.  The  factor  r  represents  a  loss 
factor,  applied  to  the  SNR,  and  caused  by  the  necessity  of 
estimating  the  noise  covariance  matrix.  The  PD  of  the  CFA.R 
detector  can  be  expressed  in  a  particularly  convenient  way  as  a 
finite  sumf 4 ' : 


p.  =  1  -TI  E  (k)('o  -  OXff  ) . 

*o  k=1 


where  L  =  K+l-N.  The  function  G  which  enters  here  is  itself  a 
finite  sum: 
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k-1  n 

V*  =  <■“'  E  £ 

n=0  ni 

In  order  to  complete  our  computat  on  of  the  PD  of  the 
likelihood  ratio  test,  we  must  take  the  expectation  value  of  this 
conditional  PD  over  the  joint  PDF  of  the  B-vectors.  These, 
however,  enter  the  final  result  only  through  the  loss  factor,  r, 
which  acts  as  a  fluctuation  model  for  the  signal.  Unlike  more 
familiar  fluctuation  models,  this  one  is  characterized  by  a 
factor  lying  in  the  range  zero  to  one.  The  present  situation  is 
similar  to  that  discussed  in  the  RMB  paper,  except  that  besides  a 
SNR  loss,  our  test  will  suffer  a  CFA.R  loss  as  well,  when  compared 
to  a  colored  noise  matched  filter  test  in  which  everything  is 
known  concerning  the  noise  or  interference. 

Although  the  loss  factor  found  here  depends  on  the  primary 
data,  through  its  B-component,  while  the  RMB  ’oss  factor  is  a 
function  only  of  the  secondary  data,  it  turns  out  that  the  two 
factors  have  exactly  the  same  PDF.  The  proof  of  this  interesting 
result  is  deferred  to  Appendix  A,  in  which  the  RMB  loss  factor  is 
also  expressed  in  our  notation,  and  the  evaluation  of  the  PDF's 
of  both  these  quantities  is  carried  out  in  parallel. 

The  PDF  shared  by  these  loss  functions  is  the  beta 
distribution: 


J(r)  =  (N+L~1)1  L 

^  *  LI  (N— 2)1  ^  ' 
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and  the  final  expression  for  the  PD  of  our  test  can  be  written 


In  this  formula/  the  H-functions  are  the  expected  values  of  the 
G's: 


1 

Hb(y)  =  J  Gk(ry)f(r)dr  . 
o 

These  integrals  are  elementary,  although  not  simple,  and  their 
detailed  evaluation  is  presented  in  Appendix  C,  where  the  result 
is  also  recast  in  a  form  more  suitable  for  computation. 
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(6) 


NUMERICAL  RESULTS  AND  DISCUSSION 


The  performance  of  the  likelihood  ratio  test  depends  only 
on  the  dimensional  integers,  N  and  K,  and  the  SNR  parameter,  a. 
The  latter  is  a  function  of  the  true  signal  strength  and  the 
intensity  and  character  of  the  actual  noise  and  interference.  Our 
analysis  deals  with  a  very  general  problem,  and  nothing  can  be 
said  about  the  anticipated  values  of  a.  The  ability  of  a  system 
to  function  effectively  in  interference  depends  principally  on 
the  arrangements  which  have  been  made  in  its  design  to  achieve  a 
good  colored  noise  matched  filter  SNR  in  its  intended 
environment.  These  arrangements  will  usually  take  the  form  of 
diversity  of  RF  inputs  in  one  form  or  another.  An  additional 
requirement  is  the  need  to  have  inputs  available  from  which  the 
actual  noise  characteristics  can  be  estimated,  and  this  is  the 
aspect  of  the  problem  which  has  been  addressed  here.  In 
particular,  for  given  values  of  PD  and  PFA,  we  can  determine  what 
SNR  is  actually  required  to  achieve  those  values  using  the 
likelihood  ratio  detector,  and  compare  that  number  to  the  SNR 
which  would  be  adequate  to  achieve  identical  performance  if  the 
noise  covariance  matrix  were  known  in  advance.  The  difference  is 
the  penalty  for  having  to  estimate  the  noise  covariance,  and  we 
expect  that  penalty  to  vary  sharply  with  the  number,  K,  of 
available  secondary  input  vectors. 

This  penalty  has  two  components:  one  due  to  the  CFAR 
character  of  the  decision  rule  and  another  due  to  the  effective 
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SNR  loss  factor.  The  latter  is  expected  to  behave  much  as  the 
results  of  the  RMB  analysis  would  predict,  based  on  the 
statistical  properties  of  the  loss  factor  alone.  The  CFAR  loss 
will  decrease  as  the  value  of  K  increases,  and  it  may  be  expected 
to  depend  largely  on  that  parameter,  while  the  SNR  loss  effect 
depends  roughly  on  the  ratio  of  K  to  N. 

These  expectations  are  borne  out  by  the  numerical 
consequences  of  our  analysis,  as  shown  in  the  accompanying 
figures.  In  Figs.  1  through  4,  probability  of  detection  is  shown 
as  a  function  of  a,  the  SNR,  for  three  detectors  (PFA  is  fixed  at 
10~6  for  these  curves).  The  detector  performing  best  is  a 
matched  filter  with  known  noise  covariance,  and  the  worst  is  the 
likelihood  ratio  detector  which,  of  course,  is  estimating  the 
noise  covariance.  The  middle  curve  (dashed)  in  these  plots  shows 
the  performance  of  a  simple,  scalar  CFAR  detector  using  L=K+1-N 
noise  samples,  and  it  differs  from  the  behavior  of  the  likelihood 
ratio  detector  only  in  that  the  SNR  loss  factor  has  been 
ignored.  This  detector  is  included  in  the  comparison  in  order  to 
show  how  much  of  the  degradation  imposed  by  noise  estimation  is 
due  to  each  of  the  two  contributing  effects. 

We  note  that  doubling  both  K  and  N  has  little  effect  on 
the  portion  of  the  degradation  due  to  SNR  loss,  while  the  CFAR 
part  is  reduced,  simply  because  K  is  being  increased.  The  curves 
also  show  the  significant  improvement  which  results  from 
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Y  of  DETECTION  vs  SNR 


PROBABILITY  of  DETECTION  vs  SNR 


PROBABILITY  of  DETECTION  vs  SNR 


increasing  the  ratio  of  K  to  N.  When  this  ratio  is  equal  to  five, 
the  SNR  loss  contribution  to  the  performance  degradation  is  about 
0.9  dB,  in  agreement  with  the  mean  value  of  the  SNR  loss  as 
obtained  from  the  beta  distribution.  Likewise,  the  CFAR 
contributions  are  directly  comparable  with  the  ordinary  CFAR  loss 
for  a  detector  of  nonf luct.uating  targets  with  no  noncoherent 
integration  (i.e,  a  single  radar  hit). 

The  detector  performance  is  characterized  in  a  different 

way  in  Figs.  5  through  8,  which  show  the  additional  SNR  required, 

when  estimating  the  noise  covariance,  to  achieve  the  same  PI)  and 

PFA  as  a  matched  filter  for  known  noise.  In  all  these  figures, 

the  Pb  is  specified  at  0.9,  but  the  results  will  not  depend 

strongly  on  tnc  chosen  PD  level,  since  the  curves  of  PD  vs  SNR 

arc  nearly  parallel  tor  the  two  detectors.  Three  Pfa  values  are 

represented  on  each  plot.  The  independent  variable  for  these 

* 

curves  is  the  number  of  secondary  vectors,  and  this  variable 
always  covers  the  range  2N  through  5N,  for  four  different  M 
values.  It  can  be  seen  that  SNR  loss  is  not  strictly  a  function 
of  the  ratio  K/N,  but  generally  decreases  with  increasing  N,  with 
this  ratio  hold  constant.  The  loss  shown  is  the  total  loss,  due 
to  the  CFAR  effect  and  the  SNR  loss  factor  itself. 

It  was  noted  earlier  that  K,  the  number  of  secondary 
vectors,  must  exceed  N,  the  dimension  of  each  of  the  data 
vectors,  in  order  to  have  a  non-singular  sample  covariance 
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matrix.  It  Is  clear  from  the  results  just  discussed  that  K  must 
exceed  N  by  a  significant  factor  if  noise  estimation  is  not  to 
cause  a  serious  loss  in  performance.  Since  N  is  the  dimension  of 
the  total  vector  of  data  used  for  detection,  the  requirements  on 
the  number  of  secondaries  can  become  very  large.  In  the  radar 
example  mentioned  in  Section  2,  N  was  the  product  of  the  number 
of  RF  channels  (M)  and  the  number  of  pulses  (J)  in  a  coherent 
processing  interval,  usually  a  large  number.  On  the  other  hand, 
our  results  are  equally  valid  for  the  case  J=l,  which  represents 
a  situation  in  which  detection  is  based  on  the  inputs  from  a 
single  pulse.  In  the  latter  case,  the  RF  inputs  could  be  the 
elements  of  an  adaptive  array,  and  N  might  then  be  a  much  smaller 
quantity. 

In  the  original  application,  the  reason  that  so  many 
secondaries  are  required  is  the  generality  of  our  formulation,  in 
which  any  interference  covariance  matrix  is  allowed.  In  the  radar 
example,  this  includes  the  possibility  of  arbitrary  correlation 
between  interference  inputs  from  pulse  returns  widely  spaced  in 
time,  although  it  is  more  realistic  to  assume  independence  (but 
not  statistical  identity)  of  the  interference  inputs  accompanying 
distinct  pulse  returns.  The  fact  that  we  have  allowed  correlation 
between  separate  pulse  returns,  but  still  assumed  that  the 
secondary  data  is  independent  of  Iho  primaries  is  somewhat 
inconsistent,  since  the  secondaries  are 


taken  irom  adjacent  range 


gate  samples.  Instead,  the  present  analysis  should  be  viewed  as  a 
formulation  of  the  problem  which  does  not  preclude  pulse  to  pulse 
independence,  but  does  not  incorporate  it  as  a  feature  of  the 
model . 

If  the  problem  is  reformulated  with  pulse  to  pulse 
independence  as  a  specific  assumption,  there  will  be  only  JM2 
real  unknown  parameters  for  the  noise  matrix  instead  of  (MJ)2 
parameters,  v/hich  could  be  an  enormous  difference.  One  expects 
the  possibility  of  improved  performance  when  fewer  unknown,  or 
nuisance,  parameters  are  being  estimated,  and  one  would  also 
expect  that  fewer  secondaries  would  be  required.  This  case  is  the 
subject  of  Part  II  of  this  study,  where  it  is  shown  that  these 
expectations  are  indeed  borne  out. 

Returning  to  the  general  problem  addressed  in  the  present 
analysis,  it  should  be  mentioned  that  the  likelihood  ratio 
decision  statistic  can  be  reformulated  in  a  way  which 
significantly  reduces  the  number  of  matrix  inversions  (or 
Cholesky  factorizations)  required.  Suppose  that  a  set  of  K+l  data 
vectors  is  specified,  and  one  of  these  is  singled  out  as  a 
primary  vector  for  the  likelihood  ratio  test.  Since  targets  can 
usually  appear  in  any  range  gate  output,  one  would  next  return 
the  selected  vector  to  the  data  pool  and  choose  another  as 
primary,  and  so  on.  Each  stage  requires  the  inversion  of  a 
different  NxN  matrix,  or  K+l  inversions  to  test  for  signals  in 
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all  the  vectors  of  the  original  set.  It  is  shown  in  Appendix  D 
that  one  can  form  a  single  sample  covariance  matrix,  using  all 
K+l  vectors,  then  invert  it  (or  factor  it),  and  use  this  matrix 
in  all  the  K+l  tests  for  target  presence  in  the  individual  range 
gates . 

Although  the  reformulated  test  is  precisely  equivalent  t  j 
the  original,  the  question  naturally  arises  as  to  the  effect  of 
signal  presence  in  more  than  one  of  the  range  gate  outputs.  This 
is  a  problem  in  any  CFAR,  and  it  is  usually  minimized  by  the 
application  of  some  screening  procedure  to  keep  signal-bearing 
vectors  out  of  the  set  used  for  noise  estimation.  Such  a 
procedure  would  limit  the  applicability  of  the  reformulated  test 
except  among  the  vectors  which  passed  the  screening  test 
themselves.  Without  screening,  the  presence  of  unwanted  signals 
in  the  secondaries  will  degrade  performance.  This  degradaton 
could  be  evaluated,  but  the  PDF  of  the  likelihood  ratio  decision 
statistic  would  become  much  more  complicated  than  it  is  without 
unwanted  signals,  and  this  topic  will  not  be  discussed  further 


APPENDIX  A 


THE  PROBABILITY  DENSITY  FUNCTION  OF  THE  LOSS  FACTOR 


The  SNR  loss  factor  derived  in  the  text  was  expressed  in 
the  form 

r  = — - — 

1  + 

where 

£t  -  fa^m  fa  • 

Before  discussing  the  PDF  of  Eg,  from  which  the  PDF  of  r 
follows  easily,  we  express  the  RMB  loss  factor,  p,  in  our 
notation.  From  the  RMB  paper, 

where  M  is  the  actual  noise  covariance  matrix,  s  is  the  signal 
vector,  and  w  is  a  weight  vector; 

w  =  kM~*s  . 

In  this  last  formula,  M  is  the  sample  covariance  matrix  of 
the  secondary  data  and  k  is  an  arbitrary  constant.  The  loss 
factor  itself  is  the  ratio  of  the  conditional  SNR  of  the  output 
of  a  filter  which  uses  w  as  a  weight,  relative  to  the  SNR  of  the 
colored  noise  matched  filter  for  known  M.  The  conditioning  in 
this  case  corresponds  to  given  values  of  the  secondary  data. 
Choosing  k=l/K,  we  obtain 

w  =  S'1  s  . 

in  terms  of  the  B  matrix  used  in  the  text,  and  then 
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(svy _ 

(*tM-1s)  (s^S~1MS_1s)  * 

Note  that  p  is  unaffected  if  s  is  changed  by  a  constant  factor. 

We  now  carry  out  the  whitening  transformation,  as  in  the 
text,  and  normalize  the  signal  vector  as  before.  The  result  is 

_(fV^)8 
_  (<*«)* 

A  unitary  transformation  is  now  applied  to  convert  the  signal 
vector  to  the  final  one  used  in  the  text,  and  the  matrices  are 
decomposed  in  the  same  way.  This  gives  the  simple  expression 

/  -  vfi 

It  is  clear  at  this  point  that  the  PDF  of  p  will  be  independent 
of  the  actual  covariance  matrix  M. 

Using  the  Frobenius  relations,  we  obtain 

(A* 


l  ^  +  '^Al'^BI  ^BA  )  » 


and  therefore 


where 


1  +  2. 


^BA 
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Note  that  the  RMB  loss  factor  depends  on  the  secondary  data  only, 
both  A-  and  B-components ,  while  r  depends  on  the  B-components  of 
both  primary  and  secondary  data. 

We  proceed  to  analyze  the  two  loss  factors  together,  and 
begin  by  conditioning  on  the  B-components  of  the  secondary  data 
vectors,  on  which  both  loss  factors  depend.  Then 

K 

?B(k)?'»(k)t 

is  a  constant  matrix,  positive  definite  and  non-singular  for  all 
sets  of  conditioning  vectors  (except  for  a  set  of  probability 
zero).  We  can  therefore  introduce  the  square  root  of  this  matrix 
and  define  the  vectors 

*  _  ..  -i/«. 
ti  “  J»  u 

and 


With  the  conditioning,  these  quantities  are  zero-mean  Gaussian 
vectors;  the  former  is  a  linear  function  of  the  elements  of  the 
B~component  of  the  whitened  primary  vector,  and  the  latter  is 
expressible  in  terms  of  the  secondary  A-cornponents ; 

t p  =  ^.rl  E  /'.(‘<)/-a0<)t  • 

k-1 

We  use  the  subscript  C  to  denote  the  present  conditioning 
and  compute  the  conditional  covariance  matrices  of  these  vectors 


1-51 


where 


Therefore 


=  As  /  /  =  ,yi* 


e.W  =  Vw*»‘  • 


K  K 

w  -  ec  E  h(kWk)f  £  ^(i)rB(Of 


E  =  #m  . 


Ec^f/  =  y,"1  . 


and  the  two  t -vectors  are  statistically  equivalent  under  the 
conditioning.  They  therefore  share  the  same  final  PDF  when  the 
conditioning  is  removed  by  averaging  over  the  secondary 
B-vectors.  Since 

and 

a,  =  (*/*,)  • 

this  proves  that  the  loss  factors  themselves  are  statistically 
identical,  and  we  continue  with  the  loss  factor,  r. 

Since  is  a  Gaussian  ( N-l ) -vector ,  the  conditional 
joint  PDF  of  its  components  is 


The  S  matrix  which  enters  here  is  itself  subject  t.o  the  Wishart. 
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PDF,  which  in  the  present  case  (in  which  the  sample  vectors  are 
of  dimension  N-l  and  the  covariance  matrix  is  equal  to  the 
identity)  takes  the  form 


In  this  formula 


f„(A) 


■  ||A||k«->i  -Tr«0 

C(N-1,K) 


C(N,K)  =  f](K-n)l 

n=1 

is  the  Wishart  normalization  factor.  The  volume  element  for  this 
PDF  will  be  written  d(A).  It  is  (N-l)2  dimensional,  ranging 
over  the  diagonal  elements  and  the  real  and  imaginary  parts  of 
tnc  upper  off-diagonals,  of  all  positive  definite  matrices,  A. 
For  our  purpose,  only  the  normalization  integral  of  the  Wishart 
PDF  is  required,  hence  we  need  not  dwell  on  the  detailed 
properties  of  this  fascinating  distribution. 

Since  depends  on  the  conditioning  data  only  through 
the  S  matrix,  its  unconditioned  PDF  can  be  written 


f .)  =  -I  t,  f -  f  !!A||  a"Tr(A* }  f *(A)  d(A) 

Tt  '  '  J  J 

As  in  the  text,  we  have  replaced  the  exponential  part  of  the 
Gaussian  PDF  by  a  trace,  this  time  involving  the  open  product 
matrix 

«  -  M.*  • 

When  we  substitute  for  the  Wishart  PDF  iri  the  expresson 
above,  we  encounter  the  integral 
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d(A)  . 


||A||K+2~N  e-Tr[A(V,  +  *)] 

This  is  the  same  as  the  normalization  integral  for  another 
Wishart  PDF,  of  dimensions  N-l  and  K+l ,  and  for  which  the 
underlying  sample  vectors  share  the  covariance  matrix 

JL  m  (IH_,  +  V)-1  . 

The  normalization  factor  for  this  slightly  more  general  Wishart 
PDF  is  just 


C(N-1,K-H)||^||k+‘ 


C(N -1,K+1)  C(N-ltK+l) 

Pm.,  +  *t'  B  [i  +  (e.fe.)]K+1 


(the  evaluation  of  the  determinant  uses  the  same  lemma  utilized  in 
Section :3)  .  Combining  these  facts,  we  ootain  the  simple  result 


/a.) = 77-— + <f, *{.))■ ^ 

1  C(N— 1,K) 


K1 


(1  +  L#)~(k+'° 


The  remainder  of  the  derivation  is  identical  to  the  final 
few  steps  given  in  the  Appendix  of  the  I<MU  paper.  The  norm  of  the 
vector  tu  is  interpreted  as  the  square  of  the  radial  coordinate 
in  a  ( 2N-2 ) -d imons i onal  Cartesian  space,  a  change  to  polar 
coordinates  is  made,  and  the  angular  coordinates  integrated  out. 
This  process  yields  the  PDF  of  and  then  a  simple  change  of 

variable  provides  the  desired  PDF  of  r: 

f(r)  - - — - (1-rf  V*'  "  - 

(K+1 -N)l(N— 2)1 
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APPENDIX  B.  EVALUATION  OF  A  DETERMINANT 

In  the  main  derivation  in  the  text,  and  again  in  Appendix 
A,  a  lemma  was  used  which  may  be  stated  as  follows: 

||A  +  0)3*11  =  ||A||(1  +  b*A-1a) 

where  A  is  assumed  to  be  nonsingular.  Because  of  this  assumption, 
we  can  write 

a  =  Ac 

and  factor  out  the  matrix  A,  and  hence  its  determinant.  It 
remains  to  be  shown  that 

l|IH  +  eb*!!  -  1  +  bfc 

where  Ijq  is  the  NxN  identity  matrix.  The  desired  evaluation 
then  follows  by  elimination  of  c. 

The  above  result,  is  proved  by  induction,  and  it  is  obvious 
for  N  =  2 .  In  general,  the  matrix 

M  3  I„  +  cb* 

is  decomposed  as  follows 

where  the  first  row  and  column  have  been  singled  out.  In  terms  of 
the  (N-l ) -vectors 


and 


c 
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we 


can  write 


w  =  b1  c 


and 


v  =  cjb 


L  =  +  cbt 


Using  the  Frobenius  relations  again,  we  have 


Of  course, 


and  wo  also  hove 


=  Mu  ||L  -  wMu"V|| 


Mu  -  1  -»•  Gfa 


wv'  =  Cjb’ob^ 


Therefore 


=  V. + 


cV 


l  1 


Assuming  the  validity  of  the  evaluation  for  the  N-l 
dimensional  case,  we  get 


=  M"(1H77‘Ti) 

\  1  H  c,  bj  / 

=  1  4  b* Cj  +  bTc  =  14  bjc 


and  the  lemma  is  proved . 
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APPENDIX  C,  COMPUTATION  OP  THE  PROP ABILITY  CF  DETECTION 


The  probability  of  detection  obtained  in  the  text  has  the 


form 


po  =  1  “TI  X  (tWo  "  l)kHk(J;) 

u 


where  L-K+l-N,  a  is  the  colored  noise  matched  filter  SNR  and  £q 

is  the  likelihood  ratio  threshold  parameter.  The  functions 

are  the  expected  values  of  the  ,  averaged  over  the  PDF  of  r, 

the  likelihood  ratio  SNR  loss  factor: 

1 


Hk(y)  =  EGk(ry)  -  J  Gk(ry)f(r  )dr  . 


The  functions  are  the  sums 


k~~i  - 

aw c 


iS)" 


(which  are  directly  related  to  the  incomplete  gamma  function), 
and  the  PDF  of  the  loss  factor  is  the  beta  distribution 

\N-2  L 


f  (r )  ~  F  (N,L)  Cl  r  /  r 


The  normalizing  factor  here  is 


F(N,l.)  =  . 

V  '  LI  (,N-2)l 


The  false  alarm  probability  is  given  by  the  very  simple 


f ormu la 
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■(if- 


and  it  can  be  seen  that  PD  reduces  to  PFA  when  a  vanishes,  since 
Gk(0)  =  l.  This  implies  that  Hj,(0)-1,  and  the  result  then 
follows  from  the  fact  that  PD  becomes  an  incomplete  binomial 


To  obtain  an  explicit  solution  for  PD,  we  substitute  for 


the  Gj<(y)  and  -write 


=  £  i„(y> 

n=0 


where 


1 

I„(y)-*|F<N.i.)  J.-"(l-r)"-i!ru"dr  . 


Next,  a  binomial  expansion  provides  the  formula 


o~ry  rL+n+pdr  _ 


The  integrals  which  enter  here  are  elementary,  and  it 
proves  useful  to  define  the  functions 


1 

T„(y)-.'  Je-'»r"dr 


In  terms  of  the  T  functions,  we  have 


I„(y>  -  £  f(N,L)  £*  (Mp2)(  -l)'Tu^(y)  , 


and  it  is  easily  verfied  that 
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The  combination  of  these  steps  provides  a  formal  solution 
to  our  problem,  but.  the  result  obtained  is  not  useful  as  a  basis 
for  numerical  evaluation.  Especially  when  y  is  small,  the  direct 
evaluation  of  the  T  functions  for  large  and  growing  m-values 
involves  the  products  of  factors,  one  of  which  is  increasingly 
large  while  the  other  is  inceasingly  small.  The  alternating  sum 
into  which  the  T  evaluations  must  be  substitued  further  stresses 
the  numerical  precision  of  the  computation. 

To  avoid  the  alternating  sum  and  its  attendant  numerical 
problems,  a  different  approach  has  been  taken  to  the  evaluation 
of  the  1  functions.  In  the  integral  which  appears  in  the 
definition  of  these  functions,  the  variable  of  integration  is 
changed  from  r  to  (l--r),  and  we  write 

I  (y)  =  ^  e-y  J_(y) 

11  ni  ” 

where 

i 

Jn(’/)  *  F(N.L)  j  ory(l-r)L+nrN-2dr  . 
o 

The  exponential  is  now  expanded  in  a  power  series  and  the 
integral  evaluated  term  by  term.  Each  of  the  integrals 


1-59 


encountered  is  the  normalizing  integral  of  another  beta 
distribution/  and  is  easily  evaluated  from  the  definition.  The  J 
functions  then  appear  as  infinite  series/  which  can  be 
expressed  in  the  form 

OO 

=  T„(y,k) 

k=0 

where 

Ki(L+n)l 
"  '  (K-tn)lLl 

and 

T  /  i  n  —  (N— 3+k)l(K+n)l  yk 
}  “(N-2)l(K+n+k)l  k!  * 

These  quantities  are  easily  generated  recursively,  and  the 
series  is  terminated  after  a  total  of,  say,  k  terms.  The 
truncation  error  is  then 


en(y-k)=DnE  Tn(y’S> 

$~k 

OO 

=  Dn  E„  Tn(y-k+s)  • 
s=0 


It  is  not  difficult  to  show  that 


■ 

£n(y>k)  =  DnTn(y-k)  E  Qn(k*S)  Ti 

S=0 

where 

Q  fk  s)  =  (N-2+k+s)l  (kln-i  k))  klsl 
n  ’  (N-2+k)l  (l'-j-n+kvs)i  (k  +  r-)l 


1-00 


In  our  application  K  exceeds  N,  and  this  ensures  that  the  0 
sequence  decreases  with  increasing  values  of  1.  Since 
Qn(k,0)=l,  we  see  that 

«n(y*k)  <  Dn  Tn(y*k)8y  * 

If  a  truncation  error  bound,  <■-,  is  prescribed,  then  the  series 
for  Jn(y)  can  be  terminated  when 

DnTn(y.k)<ce“y  . 

and  this  provides  a  simple  algorithm  for  the  computation  of  the  I 
functions.  The  sum  sequence  of  the  I's  is  the  sequence  of  H 
functions,  required  in  the  final  finite  summation  for  Pp.  The 
other  coefficients  needed  in  this  sum  (the  first  formula  of  this 
Appendix)  are  easily  generated  by  recursion.  This  procedure  has 
been  used  to  obtain  the  numerical  results  illustrated  in  the 


text. 


APPENDIX  D.  ALTERNATIVE  FORMULATION  OF  THE  DECISION  RULE 

At  first  it  seems  unlikely  that  a  test  will  perform  as 
well  if  the  primary  vector,  thought  to  contain  a  signal 
component,  is  included  with  the  secondaries  in  the  estimation  of 
the  noise  covariance  matrix.  A  simple  CFAR  example  will,  however, 
show  that  this  can  be  an  entirely  reasonable  procedure.  Suppose 
that  K+l  complex  samples  are  available,  denoted  by  z(k), 
k=0,l,,.K.  A  CFAR  test  for  signal  presence  in  z(0),  using  the 
other  samples  for  the  estimation  of  noise  level,  would  have  the 
form 

|z(0)|2>a(|z(l)f  +  -  +  |z(K)|8) 

where  the  constant  a  is  to  be  determined  by  the  assigned  value  of 
PFA.  However,  by  simply  adding 

«!*(o)|* 

to  both  sides  of  this  inequality,  and  then  dividing  through  by 
1+a,  we  find  the  equivalent  form 

|z(0)|z  (|z(0)|8  +  |z(l)f  +  -  +  |z(K)|2) 

All  the  samples  now  enter  into  the  noise  level  estimate, 
which  can  be  used  unchanged  for  tests  of  signal  presence  in  each 
of  the  other  samples  in  turn.  The  performance  attained  by  this 
detector,  and  in  particular  the  losses  caused  by  the  presence  of 
unwanted  targets  in  the  samples  used  for  noise  estimation,  is 


identical  to  that  of  the  original  form  of  the  test.  Formulas  for 
the  performance  of  such  a  CFAR  detector,  when  unwanted  targets  of 
various  kinds  are  present  in  the  noise  estimate,  may  be  found  in 
Ref.  4. 

We  expect  that  a  similar  reformulation  is  possible  in  the 
present  problem,  and  begin  by  returning  to  the  analysis  of 
Section  3.  The  likelihood  ratio  test  was  expressed  there  in  the 
form 

Min  ||  Tj  ||  0 

b 

where 

(K+1)T0  =  S  +  zz1 

and 

(K-tl)  T,  =  S  +  (z-ba)  (z-bs)f  . 

We  now  define 

nj  -f 

i>aS  +  nT 

and  note  that  this  matrix  is  K+l  times  the  sample  covariance 
matrix  formed  from  all  the  input  data  vectors. 


form 


(K+l)  Tj  =  S  -»•  (z-b8)(z-bs)t  -  ;v.zt  , 

and,  of:  course, 


(K-H)T0  =  S  . 

Next,  we  factor  out  the  new  S-rnatrix,  and  write 
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(K+l)  Tt  =  S  ^IH  +  a  b*  +  cdT^, 


where 


a  =  S  -1(z-bs) 
b  a  z— b* 


rsi  * 

csS  1  z 


and 


d  s  -z  . 

The  matrix  ratio  is  now  simply 


'i  »  __ 


I!  T0  || 


=  ||  IM  +  a  b*  +  c  d*  ||  , 


and  wc  can  again  use  the  lemma,  of  Appendix  B  in  the  following 
way.  Let  P  be  the  matrix 

P  s  IH  +  a  b*  , 


then 


||  IH  +  a  +  c  dt  ||  =  |j  P  +  c 


ii  p  ii  (i  +  (dfp  1o)j . 


Ry  the  same  lemma, 

||P|  =  l  +  (b*a), 

and  it  is  easily  verified  that 
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p"1 


a  b* 

1  +  (bTa)  ‘ 


When  these  evaluations  are  used,  we  obtain 


jp^|  =  (l+CbM^l+C^c))  -  (d*  «)(!>*  ®) 

=  (l  +  [(z“bs)f  S  -1  (z-bs)]^|l  -  (ztS-1z)j 
+  [zf  S  _1  (z-b*)]  [(z-  l>8)f  S  “l  z] 


After  this  expression  is  developed,  it  is  a  simple  matter  to 
complete  the  square,  much  as  was  done  before,  with  the  result: 

BJ(stS-lz)|8 

Q  '  ’ 

where 


LU 
II  T0|j  ' 


l  +  Q 


(•'***) 

n 


Q  *  -  (z^  S  ~l  z)^  +  |(atS~lz)|8  . 

The  minimization  over  b  is  now  trivial,  and  we  see  that  the  test 
assumes  the  form 


Q-i(sTS~lz)!8 


or 


This  is  the  desired  expression  for  the  decision  rule, 
which  is  very  much  like  its  predecessor,  but  now  involves  the 
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sample  covariance  formed  from  all  the  data  vectors.  By  using  the 
evaluation 


s-'=ts  +  z2»rl 


it  is  not  difficult, 
original  form. 


=  s's(v 

to  recover  the  li 


1  +"(?  S"fz)  / 
kelihood  ratio 


test 


in  its 
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PART  II 


(1)  INTRODUCTION 

In  Part  I  of  this  study,  a  general  problem  of  radar 
detection  was  discussed  which  was  characterized  by  the  presence 
of  unknown  non-stat ionary  interference.  The  radar  is  assumed  to 
have  a  number  of  RF  channels,  and  target  detection  is  based  on 
their  outputs  for  a  train  of  pulses  which  form  a  coherent 
processing  interval.  The  methods  of  likelihood  ratio  decision 
theory  were  applied  to  the  derivation  of  a  detection  algorithm 
for  this  problem.  Arbitrary  correlation  was  permitted  for  the 
interference  between  different  pulses,  as  well  as  between 
different  RF  channels.  A  decision  rule  was  derived  in  closed  form 
and  exact  expressions  for  the  system  performance  were  obtained. 

In  Part  1.1,  the  problem  is  reformulated  with  the  single 
additional  assumption  that  the  interference  is  statistically 
independent  from  pulse  to  pulse.  The  correlation  properties  of 
the  interference  are  still  unknown,  and  are  allowed  to  vary 
arbitrarily  from  pulse  to  pulse.  The  present  discussion  is 
largely  baseci  on  the  methods  and  results  of  Part  I,  but  the  more 
specific  problem  is  pore  difficult  to  solve,  and  we  have  not 
obtained  a  closed  form  expression  for  the  exact  likelihood  ratio 
decision  rule. 
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An  approximation  to  the  likelihood  ratio  test  statistic 
has  been  obtained  however,  and  this  will  be  derived  and  discussed 
below.  The  approximation  appears  to  be  a  reasonable  one,  and  the 
resulting  test  reduces  to  the  exact  decision  rule  for  a  single 
pulse,  and  also  to  the  likelihood  ratio  test  for  the  analogous 
multiple-pulse  problem  in  which  the  noise  covariance  matrices  are 
presumed  known.  The  form  of  this  test  itself  provides 
considerable  insight  into  the  detection  problem. 

The  exact  probability  density  function  (PDF)  of  the 
approximate  test  statistic  appears  to  be  very  difficult  to 
obtain.  An  approximation  is  developed,  however,  which  contains 
several  of  the  features  of  the  exact  PDF,  and  expressions  for  the 
probability  of  false  alarm  (PFA)  and  probability  of  detection 
(PD)  are  derived  for  this  approximation  in  the  present  report. 

An  essential  feature  of  the  analyses  of  both  Part  I  and 
the  present  Part  II  lies  in  the  assumed  existence  of  so-called 
secondary  inputs  which  can  be  used  to  estimate  the  covariance 
matrices  of  the  interference.  The  samples  from  range  gates 
adjacent  to  the  one  being  examined  for  targets  are  used  for  this 
purpose,  and  it  is  worth  specifying  in  detail  the  assumptions 
made  in  our  model  regarding  them.  It  is  assumed  that  all  external 
sources  of  interference  are  sufficiently  wide-band  in  character 
so  that,  like  internal  noise,  samples  separated  in  time  by  the 
duration  of  a  (compressed)  radar  pulse  are  independent.  The  key 
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assumption  lies  in  the  degree  of  non-stationarity  which  is 
all<  wed . 

In  the  present  model,  we  assume  that  the  covariance 
matrices  of  the  total  interference  can  be  arbitrarily  different 
for  samples  separated  in  time  by  an  interval  as  large  as  the 
pulse  repetition  interval  (PRI)  of  the  radar.  On  the  other  hand, 
the  correlation  properties  between  the  RF  channels  are  not 
supposed  to  change  so  rapidly  that  those  of  successive  samples 
(separated  by  one  pulse  length)  are  widely  different.  The 
non-stationarity  of  the  interference  is  actually  assumed  to  be 
slow  compared  to  the  pulse  length,  but  can  be  fast  relative  to 
the  PRI,  Samples  from  range  gates  adjacent  to  the  one  being 
tested  for  target  presence  can  then  be  assumed  to  be 
statistically  identical  to  the  latter,  as  long  as  they  are  not 
too  large  in  number,  while  for  successive  pulses,  the  covariance 
matrices  shared  by  this  set  of  samples  can  be  entirely 
different. 

This  model  is  admittedly  somewhat  contrived,  but  if  the 
interterence  is  allowed  to  change  arbitrarily  from  sample  to 
sample,  then  noise  covariance  estimation  is  not  possible,  and 
adaptive  detection,  in  the  sense  studied  here,  will  not  be 
feasible.  It  should  be  noted  that  the  interference  provided  by  a 
group  of  sources  with  time-varying  output  levels  does  not  result 
in  a  completely  arbitrary  time  variation  of  the  covariance 
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matrices,  and  represents  an  intermediate  case  not  addressed 
here . 

The  model  used  in  Part  I  assumed  statistically  identical 
secondaries,  but  allowed  correlation  from  pulse  to  pulse  in  both 
primary  and  secondary  inputs.  This  is  not  a  good  model  of  an 
actual  situation  which  might  be  encountered  by  a  radar,  but 
should  be  viewed  as  a  way  of  solving  the  problem  in  which  the 
pulse  to  pulse  independence  property  is  not  explicitly  employed 
in  the  derivation  of  a  decision  rule.  The  results  of  Part  I  are, 
however,  perfectly  applicable  to  the  problem  of  detection  using  a 
single  pulse. 

A  brief  summary  of  the  following  sections  of  this  report 
is  given  here.  In  Section  2,  a  likelihood  ratio  test  is  derived 
for  the  simpler  problem  in  which  the  interference  covariance 
matrices  are  assumed  to  be  known,  although  they  may  vary  from 
pulse  to  pulse.  This  represents  the  limiting  case  of  perfect 
noise  estimation,  and  the  results  provide  useful  insight  into  the 
problem  of  interest.  The  likelihood  ratio  decision  rule  for  the 
original  problem  is  derived  in  Section  3,  where  an  approximation 
is  introduced  which  is  necesary  to  obtain  a  closed  form  for  the 
test.  A  number  of  properties  of  the  approximate  test  are  derived 
in  Section  4,  and  these  tend  to  show  the  reasonableness  of  the 
approximation  made.  The  performance  of  this  test  is  analyzed  in 
Section  5,  where  a  further  approximation  is  needed  to  obtain 
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numerical  results.  The  limitations  of  the  analysis  resulting  from 
this  approximation  are  also  discussed  there.  Numerical  results 
are  included  in  Section  5,  and  a  brief  general  summary  and 
discussion  is  presented  in  Section  6.  Supporting  analysis  of  the 
final  formula  for  the  probability  of  detection  and  its  numerical 
evaluation  are  the  subject  of  the  Appendix. 


(2) 


TARGET  DETECTION  IN  KNOWN  NON -STATIONARY  INTERFERENCE 


In  the  detection  problem  addressee  in  this  study,  a  target 
is  sought  in  the  samples  corresponding  to  a  single  range  gate, 
while  K  other  signal-free  range  gate  outputs  are  available  as 
well.  If  K  is  extremely  large  and  if  all  these  signal-free 
samples  share  the  covariance  matrix  of  the  pr.mar_  sample,  we  can 
use  these  secondary  data  vectors  to  determine  an  accurate 
estimate  of  the  noise  covariance  matrices  for  each  pulse.  In  the 
limit,  these  estimates  become  the  actual  covariance's,  and  the 
problem  reduces  to  the  detection  of  a  target  in  the  primary  data, 
in  a  background  of  interference  whose  statistical  properties  are 
completely  known.  The  present  section  is  devoted  to  the  study  of 
this  simpler  problem,  which  is  useful  because  it  provides  context 
and  a  performance  comparison  for  the  more  general  problem,  and 
also  makes  a  convenient  vehicle  for  the  introduction  of  some  of 
the  notation  required. 

The  inputs  to  the  processor  are  all  data  vectors,  of 
dimension  M,  corresponding  to  the  M  RF  input  channels  of  the 
radar.  The  sample  vectors  for  a  particular  range  gate  (the  one 
being  tested  for  signal  presence)  form  a  sequence,  denoted  by 
Zj,  where  j  runs  from  1  through  J,  and  J  is  the  number  of 
pulses  in  the  pulse  train.  The  sum  cf  internal  noise  and  external 
interference  is  assumed  to  be  a  zero-mean,  circular  Gaussian 
process,  and  its  sample  vectors  for  the  pulse  have  the 
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covariance  matri> 


E0  zJ  z/  =  Mj  . 

The  symbol  Eg  stands  for  expectation  value  on  the  'noise-alone' 
hypothesis  and  the  superscript  dagger  represents  the  Hermitian 
conjugate  o£  a  vector  or  matrix.  For  different  pulses,  the  data 
vectors  are  assumed  to  be  independent,  each  pulse  being 
characterized  by  the  corresponding  covariance  matrix,  M j .  In 
the  present  secton,  all  these  covariance  matrices  are  considered 
to  be  known. 

We  use  the  subscript  '1'  to  denote  the  ‘signal-plus-noise' 
hypothesis,  and  we  characterize  the  signal  component  in  the 
following  way: 

Ei  zj  =  b  <1,  s 

where  b  is  an  unknown,  complex  scalar  amplitude  parameter,  s  is 
an  M-vector  which  represents  the  relative  signal  amplitudes  among 
the  RF  inputs,  and  the  Oj  form  a  sequence  of  complex  scalars 
which  describes  the  pulse-to-pulse  variation  in  signal  amplitude 
and  phase. 

The  signal  direction  vector,  s,  is  normalized  ao  follows 

(sT  s)  =  1  , 

i.e.  the  sum  of  the  absolute  squares  of  its  elements  is  unity.  No 
specific  normalization  convention  is  assumed  for  the  sequence 
o j .  An  example  of  the  latter  would  be  a  Doppler  progression: 
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where  6  is  the  phase  change  per  pulse  caused  by  target  motion.  It 
is  physic  aj.1^  reasonable  to  assume  that  s  is  unchanged  from  pulse 
to  pulse,  but  the  analysis  is  easily  modified  to  allow  this 
vector  to  be  a  function  of  j. 

The  joint  probaoility  density  function  (PDF)  of  the  data 
vectors  under  tne  noise-alone  hypothesis  is 


where  the  double  bar  represents 
unknown  pararneteLS  in  this  PDF. 
PDF  is  the  same,  except  that  zj 
z j  -  btfjs  in  the  exponent: 


a  determinant.  There  are  no 
Under  hypothesis  one,  the  joint 
is  replaced  by  the  quantity 


-fz.-bCTjS]1  M,  1  [zj-bffjs]^ 


ratio 


Before  maximizing  over  he 
and  take  its  logarithm: 

log  - — 


unknown 


b, 


we  form  a  likelihood 


=  M1  1  ~  t3!rbaJs]t  MJ  1  [Zj— borjS]j 

J  /  l  s 

=  ~\W  Z\°f  (sf  Mj  1  s)  +  2  Re(b*  ^  a*  (sf  Mf1  Zj)j  . 
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We  introduce  the  notation 


A,2  =  (sf  Mf1  s) 


(at  M,  z.)  (s'  M"1  z.) 

U,  s  — - — " - 2 L  > 

(s*  Mj  1  s)  Aj 

and  substitute  in  the  last  formula.  The  resalt  is 


/(b)  =  -|b|2  ^Aj2  ! of  +  2  Re(b*  £  A2  o*  Uj) 

•I  j  h  A,  a,  u.  |2 

=  —  V'  A,2  |a,!2  lb  . -  I 

Fl  i  I  A*  |<7jl2  I 

1=1  1  J 

(V  .2  *  |2 

|t,Ai  ai_ui! 

i  a/v/ 

J“i 

The  unknown  parameter,  b,  is  now  chosen  to  make  the 
quadratic  term  in  which  it  appears  vanish,  and  we  have,  simply 

!  i  Aj2  <  Uj  f 

/  -•  Max  /(b)  -  -M  . 

b  .  2  |  i2 

L  A,  ScrJ 
j-i  J 

Since  the  quantities  in  the  denominator  of  this  expression  are 
all  known,  the  likelihood  ratio  test  reduces  to 


Y'j  Aj  u*  Uj  S  >  constant. 
j=l 
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In  this  simple  problem  the  likelihood  ratio  test  is  a  'uniformly 
most  powerful  test’,  and  hence  it  yields  an  optimum  detector. 

The  quantities,  uj ,  are  recognized  as  colored  noise 
matched  filter  outputs  for  the  individual  pulses,  normalized  in  a 
particular  way.  In  the  absence  of  signal  they  have  zero  mean 
values : 

E0  uj  =  0  . 

but  when  a  signal  is  present,  we  have 

Et  Uj  =  b  tfj  . 

The  effect  of  the  denominator  in  the  definition  of  Uj  is  to 
cause  these  matched  filters  to  pass  the  signal  components  with 
unaltered  amplitudes. 

The  variance  of  the  matched  filter  output  is  easily 
computed : 

E0!Ujf  =  --4  (stMj-1  E0ZjZjt  Mj-'s)  =  ~-g  . 

Aj  3 

and  this  expression  gives  meaning  to  the  parameter,  Aj ,  as  the 
inverse  of  the  residual  noise  after  the  matched  filter  processing 
(i.e.  nulling)  on  each  pulse.  Another  way  of  viewing  the  colored 
noise  matched  filter  processing  is  based  on  a  decomposition  of 
each  data  vector  into  a  scalar  component  in  the  signal  direction 
and  an  orthogonal  component,  of  dimension  M-l.  The  orthogonal 
component  is  used  to  predict  the  noise  in  the  signal  component, 
and  this  predictor  is  subtracted  out.  The  variance  of  the  residue 
after  subtraction  is  exactly  the  inverse  of  Aj . 
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Combining  the  means  and  variances  just  obtained,  we.  see 
that  the  signal-to-noise  ratio  (SNR)  of  each  matched  filter 
output  is  just 

(SNR)j  =  |bf  Aj2  |cTj|2  , 

and  this  is  also  the  SNR  of  the  jfch  term  in  the  sum  on  which 
detection  is  based.  If  we  give  a  name  to  this  sum: 


w 


J 

E 

i=i 


-  E  V  < 


ui 


then  we  have 


E0  w  =  0 


Et  w 


f  ==  b 


J 

J=1 


E.a, 


lffi' 


and 

E0|wf  =  £  A  *  wf  E0  1 1 4  ■ 

j"i  i=1 

The  actual  SNR  of  the  likelihood  ratio  test  is  therefore  the  sum 
of  the  individual  SNR's: 

J  J 

SNR  -  |b]2  7;  A*  jor/  ^  £  (SNR)j  . 

j=l  3=1 

This  processor  is  simply  a  coherent  integrator  of  matched 
filter  outputs,  in  which  account  is  taken  of  the  varying  signal 
amplitudes  and  residual  noise  variances  of  the  individual 
pulses.  If  the  signal  amplitudes  are  constant  and  the  total  noise 
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is  stationary,  ths  processor  provides  a  gain  equal  to  J,  the 
number  of  pulses,  relative  to  the  SNR  achieved  on  each  pulse 
after  matched  filtering  (nulling).  We  will  find  analogues  of  all 
these  properties  and  parameters  in  the  more  general  problem, 
whose  analysis  begins  in  the  next  section. 
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(3)  DERIVATION  OF  THE  LIKELIHOOD  RATIO  TEST 

We  now  return  to  the  problem  described  in  Section  1  and 
assume  that  the  covariance  matrices,  Mj ,  are  all  unknown.  We 
also  introduce  the  secondary  data,  consisting  of  K  vectors  for 
each  pulse,  denoted  zj(k),  (k=l,..,K).  These  vectors  have  zero 
mean  under  both  hypotheses,  and  share  the  covariance  matrix  of 
the  primary  vector  for  the  corresponding  pulse: 

E  Zj(k)  Zj(k)*  =  M.  . 

Under  the  noise-alone  hypothesis,  the  joint  PDF  of  the  K+l 
vectors  associated  with  the  jtb  pulse  is 


fOj[zi,zj(1) . zj(K)l 


L  J.„_  e~Tr(Mf 1 T  oj) 

I  ttm  IjMjl! 


where 


T 


_1  __ 
0j  “  K+l 


zjZj* 


K 


E  zj(k)zj(k)t 

k=l 


This  representation  of  the  PDF  exactly  parallels  that  given  in 
Part  I,  and  the  corresponding  expression  in  the  signal-plus-noise 
case  is 


where 


e-Tr(Mf1Tlj(b)) 


Ttj(b)  (  (zrbajs)(zrbajs)t  +  f  ^z/k^  )  . 

k  — ■  1 

The  signal  parameters  which  enter  here  have  the  same  significance 
as  in  Section  2.  The  joint  PDF's  of  all  the  data  under  the  two 
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hypotheses  have  the  form  of  products: 


J 

0 


fo  ~  n  f°j 


and 


w-nwh)  • 

j=i 

To  form  the  likelihood  ratio  test,  each  PDF  is  maximized 
over  the  unknown  parameters  separately,  and  then  a  ratio  is 
computed.  The  unknown  signal  parameter,  b,  appears  only  in  the 
'number  one'  hypothesis,  and  we  defer  maximization  over  this 
parameter  until  last,  forming  a  ratio  first.  As  in  the  similar 
derivation  in  Part  I,  the  maximization  over  each  unknown 
covariance  matrix  is  simple;  we  obtain 

vK+l 


Max  f, 


1 


Mj  °J  (  (e7t)u  ||T0jii 


and 


Max  f.,  = 
1 


1 


\K+1 


(err)11  HT^b)!!  ) 

With  b  still  to  be  varied,  the  likelihood  ratio,  L(b),  is  a 
product  of  factors,  each  raised  to  the  power  (K+l).  It  is 
convenient  to  work  with  the  (K+l)st  root  of  the  likelihood 
ratio  instead,  and  we  write 

L(b)  »  «b)SK+1 

and 
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where 


w-n^b) . 

j=i 


^j(b)  = 


11T0,11 

llTj/b)!! 


The  sample  covariance  matrix  of  the  secondary  vectors  for 
pulse  j  is  given  by  the  sum 

K 


£  zj(k)zJ<lc)t  , 


K  k=l 

and  this  quantity  is  also  the  maximum  likelihood  estimator  of  the 
unknown  covariance  matrix  for  this  pulse,  based  on  the  secondary 
data  alone.  In  terms  of  this  matrix,  we  have 

(K+l)T0j  =  KMj  +  Zjz/  . 
and  by  the  determinant  lemma 


|| A  +  a btH  =  ||A||  |l  +bfA  'aj  , 


for  which  a  proof  is  given  in  Part  I,  we  obtain 


(K+1)MI|T01||  =  KU  liMjll  (l  +  . 


1  +  ^(z/M,-1z1) 


By  the  same  reasoning  we  find 

(K+l)M|lT„(b)||  =  K“||Mj||  (l  +  I  [Zj-bo-jsy  M,-1  [zrb<7jS]^ 

for  the  signal-plus-noise  case,  and  thus  we  have  the  ratio 

<j(b) 


1  +  ^  [Zj-bo-jS^Mj  ^Zj-bo-js] 


In  analogy  to  the  quantities  introduced  in  Section  2,  we 


make  the  definitions 


and 


Uj  (sf  Mj"1  s)  Aj2 

It  should  be  noted  that  these  new  variables  involve  the  estimator 
of  the  covariance  matrix,  and  not  the  unknown  Mj  itself.  It  is 
also  convenient  to  use  the  temporary  notation 

With  the  help  of  these  definitions  we  can  make  the 
evaluation 

[xpbcrJsjtMj~1[Zj---bCTjS]  =  Nj  -  2  Relb"  Af  t  |'bi  iffj'f  Aj 

==  Aj2  lbcrj  -  u/  +  Nj  -  Aj2|Uj|2  . 

The  likelihood  ratio  for  the  jth  pulse  is  then  given  by 


<j(b)  “ 


V  K.  N1 . . . . . 

1  +  W  ■  k xjV, - ai 


Next,  we  introduce  the  quantity 


and  eliminate  Uj,  so  that 


1  -1  >-g  H  ]}{  Aj2!u/' 


1  1  ^  Aj  \Uoi  -  Uj! 
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The  signif  .cance  of  Ij  can  be  seen  by  defining  an  inner 
product  as  follows: 

(a.b^^M^b)  . 

Here,  a  and  b  are  arbitrary  complex  M-vectors,  and  the 
corresponding  norm  is 

||  af  »(a'Mj-‘a)  . 

t 

This  definition  is  possible  because  the  sample  covariance  matrix 
and  its  inverse  are  positive  definite  with  probability  one  (so 
long  as  K  exceeds  M,  a  condition  we  now  impose).  We  can  now  write 

A*  =  (is  f 

and 


(s.Zj) 

II  si! 


This  last  expression  is  the  component  of  Zj  in  the  direction  of 
the  'unit  signal  vector' 


s 


The  portion  of  zj  which  is  orthogonal  to  the  signal 
direction  (in  the  sense  of  this  new  inner  product)  is 


and  in  terms  of  this  projection  we  obtain  the  desired  expression 
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K2j  =  Nj  —  Aj2|ujf 


Thus  KEj  is  the  norm  of  that  portion  of  the  data  vector  which 
is  orthogonal  to  the  unit  signal  vector.  This  separation  of  a 
data  vector  into  a  'signal'  component  and  an  orthogonal 
component,  relative  to  some  noise-dependent  inner  product,  will 
appear  again  in  connection  with  the  coherent  integration  of 
pulses. 


It  will  be  shown  in  a  later  section  that  Ej  is  exactly 
analogous  to  the  quantity  %  of  Part  I,  and  also  that  the 
variables 


1 

r  = - 

J  1  +  E, 

represent  SNR  loss  factors,  described  statistically  by  the  Beta 
distribution . 

Returning  to  our  likelihood  ratio,  we  divide  numerator  and 
denominator  by  (1  +  Ej)  and  make  use  of  the  above  definition  to 
obtain  the  formula 


1+  il-jA,  ib<T,-Uj 

The  final  likelihood  ratio  is  the  product  of  J  of  these  factors, 

A  /V 

maximized  over  b,  and  we  note  that  r j ,  Aj  and  uj  are 
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A  i  m  ft  dtaitni  n  Lf»‘i  H  UijAyi 


4 


sufficient  statistics  for  the  decision  process.  These  scalar 
parameters  are,  in  turn,  simple  functions  of  the  norms  of  s  and 
zj,  and  the  inner  product  (s,zj),  as  defined  above. 

We  cannot  derive  an  exact  likelihood  ratio  decision  rule, 
because  the  required  maximization  of  the  product  of  factors  (over 
b)  cannot  be  carried  out  in  closed  form.  However,  we  note  that 
the  parameter  K,  the  number  of  secondary  data  vectors,  will 
control  the  accuracy  of  the  noise  covariance  estimate,  and  we 
know  that  K  will  have  to  be  large  compared  to  M  in  order  to  avoid 

/v 

significant  SNR  loss.  For  large  K,  the  random  variables  Aj  and 

A 

uj  will  tend  toward  the  constants  Aj  and  u j ,  as  the  sample 
covariance  matrix  tends  to  the  true  noise  covariance  matrix.  In 
the  same  limit,  the  loss  factors,  r j ,  will  approach  unity 
because  of  the  properties  of  the  Beta  distribution  to  which  they 
are  subject.  Thus  all  the  quantities  in  the  likelihood  ratio, 
other  than  K  itself,  remain  bounded  as  K  increases  and  we  will 
assume  that  for  practical  values  of  K,  both  the  numerator  and 
denominator  of  the  likelihood  ratio  (for  each  pulse)  take  the 
form  of  unity  plus  a  small  quantity. 

Motivated  by  this  reasoning,  we  make  the  approximation 

J 

/(b) = n  ^b) 

M 
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products  of  a  unitary  matrix  and  a  new  set  of  random  vectors. 
These  new  random  vectors  are  statistically  indistinguishable  from 
their  predecessors,  and  it  would  only  be  confusing  to  introduce  a 
new  notation  for  them.  Tracing  this  transformation  through  the 
inner  products,  we  find  that  only  the  normalized  signal  vector  is 
changed:  t  is  replaced  by 

a  ty 

where  U]_  is  the  unitary  matrix  characterizing  this  last 
transformation.  Any  unit  vector  in  the  complex  N-space  can  be 
realized,  as  ti ,  by  such  a  transformation.  In  particular,  we 
can  cause  ti  to  be  a  'coordinate  vector1,  for  which  a  single 
element  is  unity,  the  remaining  (N-l)  elements  vanishing.  It  is 
for  this  reason  that  the  PDF  of  n  depends  on  M  only  through  the 
meaning  of  the  signal  amplitude  parameter,  b.  In  fact,  this  PDF 
can  depend  only  on  b,  N  and  K,  and  hence  the  false  alarm 
probability  of  the  like] ihood  ratio  detector,  namely 

PFA  =  Prob{77  >  7]0  \  , 

is  independent  of  M,  and  this  is  the  generalized  CFAR  property 
claimed  in  section  1. 
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(4)  PROPERTIES  OF  THE  APPROXIMATE  LIKELIHOOD  RATIO  TEST 

The  approximate  likelihood  ratio  (ALR)  test,  derived  in 
the  previous  section,  has  a  number  of  properties  which  suggest 
that  it  is  a  reasonable  detection  procedure,  in  spite  of  the 
approximation  made  in  its  derivation.  Before  discussing  these,  it 
is  perhaps  worthwhile  to  point  out  just  what  is  involved  in  its 
actual  implementation.  The  procedure  requires  the  estimation  of  a 
covariance  matrix  for  each  pulse,  carried  out  by  averaging  the 
open  (dyadic)  products  of  the  secondary  data  vectors.  Two  sets  of 
linear  equations  must  then  be  solved,  which  we  may  write  in  the 
form 


Mjqj  =  Zj  . 

The  three  inner  products  are  then  given  by 

(s^j^s)  =  (vr^s) 

(stMJ~1zJ)  =  (w/zj) 

and 

(z/  Mj-1  Zj)  =  (q Zj)  • 

The  variables  rj ,  Aj  and  uj  are,  as  we  have  seen,  defined 
directly  in  terms  of  these  inner  products,  and  only  these 
variables  are  needed  for  the  construction  of  the  ALR  test. 


11-22 


By  far  the  greatest  task  is  the  sclution  of  the  linear 
equations,  normally  carried  out  by  means  of  a  Cholesky 
factorization  technique.  A  nulling  scheme  which  employs  the 
colored-noise  matched  filter  method  directly,  and  which  is 
designed  to  compute  new  weights  for  each  pulse,  will  have  to 
solve  the  first  of  the  above  equations  to  find  the  appropriate 
weight  vector.  After  factorization  of  the  sample  covariance 
matrix,  it  does  not  require  much  extra  computation  to  solve  the 
second  equation  as  well.  The  inner  products  are,  of  course, 
‘simple  complex  sums,  as  are  the  sums  over  j  which  enter  into  the 
test  itself.  Therefore,  the  implementation  of  the  ALR  decision 
rule  would  require  only  a  moderate  increase  in  computational 
complexity  over  such  a  matched-filter  nulling  processor. 

From  the  form  of  the  ALR  test,  and  the  definitions  of 
rj,  Aj  and  uj  ,  it  can  be  seen  that  the  test  is  unaltered  if 
the  signal  vector,  s,  is  changed  by  a  scale  factor.  Thus  the 
detection  rule  is  not  affected  by  our  choice  of  signal  vector 
normalization.  More  obviously,  the  test  is  unchanged  if  the 
sequence  Oj  is  multiplied  by  a  common  scale  factor  as  well. 
Finally,  the  test  is  invariant  to  a  scaling  of  all  the  data 
vectors,  primary  and  secondary,  by  a  common  scalar  factor.  This 
last  represents  a  weak  form  of  CFAR  behavior.  It  will  be  shown 
later  that  the  test  is  not  a  true  CFAR,  in  the  sense  that  its 
probability  of  false  alarm  (PFA)  is  completely  insensitive  to  the 
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actual  covariance  matrices  of  the  noise  on  all  of  the  pulses. 
However,  it  will  turn  out  that  the  test  is  approximately  a  true 
CFAR,  and  we  will  be  able  to  get  some  idea  of  the  nature  of  its 
departure  from  the  desirable  true  CFAR  performance. 

If  K,  the  number  of  secondaries,  becomes  very  large,  the 
ALR  test  passes  over  into  the  likelihood  ratio  test  derived  in 
Section  2,  for  the  case  of  known  noise.  This  follows  from  the 
convergence  of  the  noise  covariance  estimators  to  the  true 
covariance  matrix  in  this  limit,  and,  as  discussed  before,  the 
convergence  of  the  Beta-distributed  loss  factors,  r j ,  to 
unity.  In  the  final  form  of  the  ALR  test,  given  as  the  last 
equation  of  Section  3,  the  summation  over  j  on  the  right  side 
will  remain  finite  as  K  increases,  and  hence  this  sum,  divided  by 
K,  will  tend  to  zero.  If  the  quantity 

U-\ 

k  ~r 

is  then  defined  as  a  new  threshold  constant,  the  limiting  form  is 
seen  to  be  identical  with  the  known-noise  test  of  Section  2. 

In  the  case  of  a  single  pulse  the  ALR  test  is  exact,  since 
the  approximation  we  have  made  is  irrelevant  in  this  situation. 
This  case  is  actually  the  same  as  the  general  problem  treated  in 
Part  I,  since  we  have  made  use  of  no  special  structure  for  our 
M-vectors.  If  we  put  J  =  1,  the  ALR  test  becomes 

rA2jGf  >K^1  +  |rAa|Sf) 
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or 


1  0  2 1^|2  .  ,  . 
rA  ju]  >  -  1  . 


K 

We  have  dropped  the  subscript  j  in  these  expressions. 

A  A 

Substituting  for  r,  A  and  u,  we  obtain 

|(sf  M_1  z)|2 

>  (*o  -  D(1  +  E)  . 

K  (s' M  ‘s) 

again  without  the  subscripts.  Finally,  substituting  for 
find 

K(stM  *s)  Zo  \  K  > 


we 


which  is  identical  to  the  decision  rule  derived  in  Part  I. 

This  test  can  be  put  in  another  form  which  will  provide  an 
interesting  analogy  to  the  multiple-pulse  ALR  decision  rule. 

Using  the  inner  product  notation  introduced  in  Section  3  and  the 
definition  of  the  orthogonal  component  of  the  data  vector  with 
respect  to  that  inner  product,  we  can  write  the  next  to  last  form 
for  the  single-pulse  test  as 

'fs  z)i2  /  \ 

L  * --->(<0-1)(K  +  (zllZl))  • 

(s,s)  \  / 

Returning  to  the  multiple-pulse  test,  we  define  a  new  inner 
product,  using  square  brackets  to  distinguish  it,  as  follows: 

J 

[a-b]  =  rjAj2aj*bj  . 

j=l 

Here  a  and  b  are  arbitrary  sequences  of  complex  numbers,  or 
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J-vectors.  Using  this  inner  product/  the  ALR  test  can  be  written 

|[cr,u]|2  /o-lLr-i 
[a, a]  \ 

Again  we  can  separate  the  data  vector/  this  time  the 

A 

J-vector  u,  into  two  components  relative  to  the  signal  vector, 
a,  defining  an  orthogonal  component  as  follows: 

-  - 

u  BU-  C  . 

[a.crj 

/v 

The  norm  of  u  is  then  the  sum 

o  T  |2 


r_,  H>.u]r  ~  , 

[u,u]  =  - - -  -  +  [uruA] 

[CT.crJ 


and  the  ALR  test  reduces  to 


~  >  U0  ~  MK  +  RivS)  ’ 

,a\  \  1 


!l>.u]|' 
[> 


which  stands  in  remarkable  analogy  to  the  single-pulse  form  of 
the  test,  just  derived  above. 

To  procede  with  the  analysis  of  the  ALR  test  we  follow  the 
procedure  of  Part  I  and  introduce  a  change  of  variables,  to 
'whitened'  coordinates.  Suppose  the  actual  noise  covariance 
matrix  for  the  data  vectors  of  the  jt-*1  pulse  is  M j .  New 
vectors  are  defined  by  the  equations 


s  Mf2zj 


and 


7j(k)3Mj  2Zj(k)  . 
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The  covariance  matrix  shared  by  all  K+l  cf  these  vectors  is  the 
M *M  identity  matrix,  and  the  same  is  true  for  the  data  vectors 
for  all  J  pulses.  The  sample  covariance  matrix  of  the  whitened 
secondaries  is  given  by 

K  l  l 

UCj-ty;  ?j(k)?i(k)*  -  Mj-SMjMj-8  . 
k=l 

The  signal  vector  is  also  changed  by  the  whitening 
transformation  into  a  different  vector  for  each  pulse: 

tj  s  Mf2S  , 

and  here  we  see  it  would  be  very  simple  to  allow  the  original  s 
vector  to  be  a  function  of  j.  The  ordinary  norm  of  the 
transformed  signal  vector  is 

(t/tj)  =  (sfMj"ls)  =  Aj2  . 

the  quantity  encountered  in  Section  2,  which  is  the  inverse  of 
the  residual  noise  that  would  remain  after  nulling  with  a  filter 
matched  to  the  actual  noise. 

In  Part  I  we  were  able  to  choose  a  convenient  signal  norm 
which  made  this  quantity  unity,  since  there  was,  in  effect,  only 
one  pulse.  That  choice  cannot  be  made  here,  and  hence  a  simple 
norm  has  been  chosen  for  s  itself.  The  presence  of  the  Aj  in 
our  expressions  will  be  the  chief  difference  between  the  present 
analysis  and  that  of  Part  I.  The  evaluations  that  are  made  next 
follow  very  closely  those  of  Part  I,  and  only  enough  detail  will 
be  given  to  make  clear  the  small  differences.  Besides  the 
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appearance  of  th i  Aj ,  we  are  now  working  with  the  sample 
covariance  matrices  themselves,  instead  of  the  S-matrices  of  Part 
I,  which  differ  by  a  factor  of  K. 

We  introduce  the  inverses  of  the  whitened  sample 
covariance  matrices: 

?.  -  Uif1  =  M,_1  Mj  2 

and  substitute  the  definitions  in  our  three  inner  products.  We 
obtain 

and 

(Zj^  Mj  1zi)  =  • 

The  new  random  vectors  are,  of  course,  still  Gaussian,  and 
the  secondaries  still  have  zero  means  under  either  hypothesis. 

The  whitened  primary  vector  has  zero  mean  in  the  absence  of 
signal,  and  mean  value 

E,,,  =  MjhE12j  =  b<ritj 

under  the  1 s ignal-plus-noi se 1  hypothesis.  These  variables  are 
thus  completely  characterized  statistically,  and  we  note  that  the 
actual  noise  covariance  matrices,  Mj ,  appear  only  in  the  new 
signal  vectors,  the  t-j. 

Following  the  method  of  Part  I,  we  next  make  unitary 
transformations,  one  for  each  pulse,  which  will  leave  all  the 
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Gaussian  random  vectors  statistically  unaltered  except  for  the 
mean  values  of  the  primary  vectors.  These  mean  values  are  • 
proportional  to  the  transformed  signal  vectors,  and  we  choose  the 
unitary  transformations  to  make  each  signal  vector  proportional 
to  the  'coordinate  unit  vector'  whose  transpose  is  [1,0,...  ,0], 
Since  the  norms  of  the  signal  vectors  are  unchanged  in  this 
process,  the  transformed  tj  will  be  given  by 


In  this  vector,  the  'one'  is  a  scalar  and  the  'O'  is  an  (M-l) 
vector.  We  can  think  of  these  unitary  tranf ormations  as  being 
combined  with  the  preceding  whitening  transformations  (which  were 
• n  any  case  not  unique),  and  hence  we  make  no  change  in  the 
notation  to  reflect  the  unitary  transformations. 

We  see  now  that  the  matrices.,  M j ,  affect  the  statistical 
character  of  the  inner  products  only  through  the  numbers,  Aj ,  a 
single  scalar  for  each  pulse.  Moreover,  these  inner  products  are 
sufficient  statistics  for  decision,  since  the  ALR  test  itself 
depends  on  the  data  vectors  only  through  them.  The  only  other 
parameters  which  enter  the  test  are  the  signal  amplitude 
constants,  the  o j ,  and  we  shall  find  that  the  test  statistics 
depend  only  on  the  products,  Aj Oj . 

The  data  vectors  and  sample  covariance  matrices  are  now 
decomposed  into  blocks,  separating  out  the  'one'  component. 
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L  now  depends  on  the  dimension,  M,  of  the  data  vectors  for  each 
pulse,  instead  of  the  total  number,  MJ,  of  data  vectors  which 
enter  in  the  detection  process.  This  difference  is  due,  of 
course,  to  the  more  specific  assumption  made  here  concerning  the 
covariance  structure  of  the  interference,  namely  pulse-to-pulse 
independence , 

When  all  these  substitutions  are  carried  out,  new 
combinations  of  variables  appear,  which  motivate  the  definitions 

yj  s  fjA  ~  ^jB 


y3(k)  =  ?jA(k)  ~  *^JAB^jBB  J'jB(k)  • 

These  variables  are  the  precise  analogues  of  y  and  y(k), 
encountered  in  Part  I,  and  we  also  obtain  the  relation 

K 

*^jAA  “  ^jAB-^jBB  -^jBA  “  ft  ^  ’ 

Note  the  appearance  of  the  factor  1/K  in  this  formula,  which 
results  because  the  left  side  is  expressed  in  terms  of  the  sample 
covariance  matrices  themselves. 

When  conditioned  on  the  B  components  of  the  primary  and 
secondary  vectors  of  pulse  j,  the  yj  and  yj(k)  are  Gaussian, 
and  yj  is  independent  of  the  yj(k).  The  K*K  conditional 
covariance  matrix  of  the  yj(k)  has  the  same  structure  found 
earlier;  it  is  an  idempotent  matrix  with  L  unit  eigenvalues,  the 
remaining  M-l  eigenvalues  being  zero.  Conditionally,  the  means  of 
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the  y  j  ( k )  are  always  zero,  and  that  of  y  -j  is  given  in  general 
by 


The  subscript  B  denotes  the  conditioning,  and  obviously  this  last 
mean  is  zero  in  the  absence  of  signal.  Finally,  the  conditional 
variance  of  yj  is  related  to  the  loss  factor,  as  follows: 

EB|yJ-bAJffJ!8  =  l  +  E,=^  • 

In  terms  of  these  y-variables,  it  is  easily  shown  that 


u3  = 


_yj 


and 


2  2  _  K  Ai 

A  • 


y  |v  f\r\f 

k=l 


The  variables  which  enter  into  the  three  basic  inner  products 
have  now  been  written  in  terms  of  the  y's  and  E j ,  whose 
statistical  properties  have  been  completely  described.  The 
'square-bracket'  inner  products  themselves  are  expressible  in 
terms  of  these  quantities,  and  we  find 

J  J 


3=1  3=1  L  ly.(k)l 

k°l 


[a.a] 


rjA]2kji2 

£  |yj(k)!2 

k=>l 
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and 


b.5]  =  £  r.WVKE 

Fi  J"i 


r|Vi*y  1 

5  |y,00f 

k=l 


We  have  already  seen  that  the  ALR  test  takes  a  very  simple  form 
when  written  in  terms  of  these  inner  products. 

Following  the  analysis  of  Part  1,  we  find  that  the  sums 


K 

Tis  Z  iyj(k)!2 

k=l 

are  Chi-squared  random  variables,  independent  of  conditioning  on 
the  B-components ,  and  each  Tj  is  subject  to  the  PDF 


(L~ 1)! 

Note  that  the  number  of  degrees  of  freedom  is  related  to  L  and 
not  K,  because  of  the  mutual  correlation  of  the  yj(k).  We  now 
introduce  new,  normalized  random  variables  instead  of  the  T  j  , 
namely 


T 


3 


These  new  variables  have  unit  mean  values,  and  each  has  variance 
equal  to  1/L.  These  properties  are  central  to  an  approximation  we 
shall  make  in  the  next  section,  in  connection  with  the  performance 
of  the  ALR  test.  Like  all  variables  relating  to  different  pulses, 
the  Tj  are  mutually  independent. 
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Again  following  Part  I,  we  replace  the  yj  by  new 
variables,  '  j,  which  are  defined  so  that  they  have  unit 
variance,  conditioned  on  the  B-components : 


(1  + 


-  —  r,%' 


The  loss  factors  now  appear  in  the  conditional  means: 

EBwj  =  brjSAJal  • 

To  express  the  ’square'  inner  products  in  compact  form,  it 
is  natural  to  use  the  definition 

s  rjZ  AjaJ  • 

These  quantities  are  random  variables,  but  they  are  dependent 
only  on  the  B-components  of  the  data  vectors.  This  last 
definition  allows  us  to  write 

EBwJ  ^  b^i 

and  the  inner  products  are  then  given  by  the  equations 

1  p  '•i  1  y  ^ 

K  [“•»]  "  L  £  — 

-  Ter  <rl  -  -  V 
K  L  pi  rj 


1  r  -T  1  V  ^ 


Wj 


The  ALR  test  itself  is  therefore  statistically  equivalent  to  the 


test 
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A  I/*, f 

]=1  Tj 


1  +  f  £  T 

L  j„i  Tj 


Of  course,  the  transformations  we  have  made  were  based  on  the 
actual  noise  covariance  matrix,  and  the  form  of  the  test  just 
given  will  be  used  only  to  study  its  performance;  it  is  not  an 
alternative  representation  to  the  original  form,  which  was 
written  in  terms  of  the  observables  themselves. 

Conditioned  on  the  B-coraponents  of  the  whitened  vectors, 
the  uj  are  constants  and  the  wj  are  simple  Gaussian 
variables.  We  note  that  the  conditioning  is  represented  now  only 
by  the  presence  of  the  loss  factors,  hence  we  can  interpret  the 
conditional  probabilities  as  being  conditioned  on  a  set  of  values 
of  these  J  independent  random  variables.  Besides  these,  only  the 
products  Aj  a j  appear  in  the  test,  and  also  in  the 
statistics,  since  the  conditional  means  of  the  wj  are  equal  to 
the  numbers  u j,  which  are  proportional  to  these  same  products. 

Because  of  the  presence  of  the  products  Aj  aj  in  the 
decision  rule,  it  is  not  a  true  CFAR,  in  the  sense  that  the 
probability  of  false  alarm  is  not  totally  insensitive  to  the 
actual  noise  covariance  matrices.  However,  the  PFA  is  invariant 
to  any  permutation  of  these  numbers,  since  the  random  variables 


which  enter  the  test  are  independent  from  pulse  to  pulse  and 
statistically  indentical  for  different  pulses.  The  PFA  is  also 
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unchanged  if  all  the  Aj  are  changed  by  a  common  factor,,  We 
conclude  that  the  PFA  depends  on  the  variability  of  the  sequence 
of  products,  Ajtfj,  in  some  normalized  way,  such  as  the  ratio 
of  standard  deviation  to  mean  of  these  numbers.  It  seems  likely 
that  the  dependence  of  the  PFA  on  this  variability  is  not  strong, 
and  in  the  next  section  an  approximation  will  be  introduced  which 
causes  this  dependence  to  disappear  entirely.  This  approximation 
can  be  used  to  set  a  unique  threshold  for  the  test,  for  a  given 
assigned  PFA,  and  the  actual  PFA  will  then  vary  somewhat  from 
this  assigned  value,  according  to  the  actual  pulse  to  pulse 
variations  in  the  level  of  the  residual  noise  after  nulling. 


(5) 


PERFORMANCE  OF  THF  APPROXIMATE  LIKELIHOOD  RATIO  TEST 


To  evaluate  the  performance  of  the  ALR  test  exactly,  it  would  be 
necessary  to  obtain  the  PDF  of  the  quantity 


The  probability  that  the  detection  threshold  is  exceeded  is  then 
equal  to 

Prob(£>0)  . 

Evaluation  of  this  probability  under  the  noise-alone  hypothesis 
yields  the  PFA  of  the  decision  rule,  and  under  the  signal-plus- 
noise  hypothesis  we  obtain  the  PD.  We  have  shown  that  the  ij 
are  independent,  normalized  Chi-squared  variables’,  independent  of 
any  conditioning,  and  that  the  wj  are  independent.  The  wj  are 
also  Gaussian,  when  conditioned  on  the  loss  factors  contained  in 
the  uj .  Finally,  the  loss  factors  themselves  are  independent 
and  satisfy  the  Beta  distribution  given  in  Section  4.  We  would 
like  to  be  able  to  evaluate  the  conditional  PDF  of  K,  and  then 
remove  the  conditioning  by  taking  the  expectation  value  of  this 
probability  with  respect  to  the  loss  factors.  This  procedure  was 
feasible  for  the  general  problem  analyzed  in  Part  I,  but  appears 
to  be  intractable  here. 

It  was  pointed  out  in  Section  4  that  the  PDF  of  the 
Tj  has  mean  value  unity  and  variance  equal  to  1/L,  where 
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L=K+1-M.  Since  L  will  have  to  be  large  compared  to  unity  in  order 
to  control  SNR  losses  (i.e.  to  keep  the  loss  factors  close  to 
unity),  it  may  be  expected  that  the  Tj  will  not  differ  yreatly 
from  unity  themselves.  A  family  of  PDF  curves  for  this  normalized 
Chi-squared  distribution  is  presented  in  Fig.  1,  for  various 
values  of  L.  It  can  be  seen  that  the  PDF  becomes  relatively 
narrow  and  quite  symmetrical  about  the  mean,  for  values  of  L  in 
excess  of  about  50.  This  suggests  making  the  simplifying 
approximation 

Tj  »  1  . 

in  the  expression  for  K,  which  then  becomes 


This  approximation  greatly  simplifies  the  PDF  of  5,  and  it 
is  only  this  simplified  form  that  will  be  discussed  here.  The 
penalty  associated  with  the  simplification  is,  of  course,  a 
failure  of  our  results  to  describe  all  aspects  of  the  original 
problem.  In  Part  I  it  was  shown  that  the  likelihood  ratio  test 
derived  there  was  equivalent  to  a  simple  CFAR  detection  process, 
in  which  the  threshold  was  estimated  from  secondary  data,  and  the 
target  exhibited  a  Beta-distributed  fluctuation.  The  threshold 
estimation  in  that  problem  was  represented  by  a  single  variable 
exactly  like  one  of  our  t j .  The  resulting  performance 
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analysis  (in  Pari.  I)  also  showed  that  the  total  SNR  loss,  could  be 
recognized  as  the  sum  of  two  contributions:  one  due  to  the  target 
fluctuation  and  another  due  to  the  threshold  estimation.  The 
latter  effect,  which  provided  the  CFAR  behavior,  accounted  for  a 
SNR  loss  very  close  to  the  standard  CFAR  loss  for  a  non¬ 
fluctuating  target  and  threshold  estimation  by  summation  of  noise 
sample  powers.  The  number  of  terms  entering  in  this  summation  is 
exactly  L.  In  the  present  problem  there  is  a  Tj  factor  for  each 
term  in  the  sums  which  enter  in  the  decision  statistic,  since 
noise  estimation  is  carried  out  for  each  pulse.  Each  pulse  will 
then  show  a  CFAR  loss,  and  we  expect  that  the  overall  performance 
of  the  actual  ALR  test  will  have  a  CFAR  loss  of  the  same  order  of 
magnitude  (in  db)  as  that  for  a  typical  pulse.  It  is  this  aspect 
of  performance  that  will  not  be  described  by  our  results  when  the 
simplifying  approximation  is  made. 

In  the  approximating  expression  for  5  we  can  interpret  the 

Wj  and  Uj  as  components  of  complex  J-vectors  and,  as  before, 

make  use  of  the  inner  product  notation.  Then  we  have 

J 

J=1 

J 

Y,  iMjf  =[/*./*] 

1=1 

and 

J 

V  kj!2  =  [w(w]  . 
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J 

{'  =  KT  -  -  ('o-DL 

and  note  that  the  probability  of  exceeding  the  detection 
threshold  is  equal  to 

Prob(£’>0)  . 
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The  components  of  w'  are  independent  and  Gaussian, 
the  loss-factor  conditioning,  and  each  has  unit  variance, 
conditional  mean  of  the  vector  w*  is 


g  iven 
The 


EBvr'  =  b/x'  , 

and  therefore  all  components  of  w’  except  the  first  have  zero 
mean.  The  latter  component  has  mean  value 

r  -t  1  1 

EBw'1  =  b[/x./x]2  =  a2 

where 

» s  n>C  = ihf*  £  riAiW  • 

J=1 

This  last  parameter  is  the  basic  signal  parameter  of  the 
simplified  problem,  and  we  see  that  it  is  just  like  the  SNR  of 
the  known-noise  test  (see  Section  2)  except  for  the  presence  of 
the  loss  factors,  the  r j . 

In  the  noise-alone  situation,  a  =  0,  and  the  simplified 
test  performance  depends  only  on  L  and  the  threshold  parameter; 
it  is  in  this  approximation  a  true  CFAR .  When  a  signal  is  added, 
its  presence  is  felt  only  through  the  parameter,  a,  and  we  shall 
calculate  the  conditional  PD  for  this  case,  as  a  function  of  a. 
It  will  not  be  possible  to  average  this  PD  over  the  PDF  of  a,  in 
order  to  remove  the  conditioning,  because  of  the  difficulty  of 
dealing  with  a  weighted  sum  of  Beta-distributed  variables.  Again 
we  appeal  to  the  results  of  Part  I,  where  it  was  shown  that  the 
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exact  SNR  loss  caused  by  the  loss  factor  was  very  similar  to  the 
expected  value  of  the  loss  factor  itself.  We  therefore  note  that 


J 

E a  =  E r |bf  £  A*|<»/ 

J=1 

where 

Er  =  J'  rf(r)dr 
o 

=  L+!  „  ?lt_2rN 
N+L  K+l 

is  the  expected  loss.  The  expected  value  of  a  is  thus  equal  to 
the  known-noise  SNR  value  times  the  expected  value  of  a  typical 
loss  factor,  as  found  from  the  Beta  distribution. 

It  remains  to  derive  the  probability  of  threshold  crossing 
for  fixed  a,  from  the  simplified  version  of  the  test  statistic. 

We  define  the  characteristic  function  of  the  random  variable 

$(A)sEeiXr  . 

and  note  that  the  desired  probability  can  be  expressed  as  an 

integral  in  the  complex  A-plane,  as  follows 

--1* 

Prob(f>0)«-t  J  *(«“ 

The  interchange  of  order  of  integration  required  to  obtain  this 
result  is  made  valid  by  the  displacement  of  the  contour  below  the 
real  axis  in  the  complex  plane.  The  characteristic  function 
itself  can  be  expressed  in  the  form 

$(A)  =  e"iUo-1)LX  * 
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where 

and 

$2(X)  =  EeiXJ4|v'j|2 

is  the  characteristic  function  of  the  sum 

s  -  £  Kjf  . 

j=2 

Since  S  is  the  sum  of  the  squares  of  (J-l)  complex 
Gaussian  variables,  each  with  mean  zero  and  variance  unity,  it  is 
a  Chi-squared  variable  of  2J-2  degrees  of  freedom,  and  its 
characteristic  function  is  therefore 

*sw  =  (»  -  • 

The  random  variable  w'j_  is  conditionally  Gaussian,-  with  unit 
variance  and  a  mean  squared  value  of  a,  hence  its  characteristic 
function  is 


Substituting , 


e 

r-Tx 


we  obtain 

a 

el-iX 

^xSTiU0-ijxpf 


$(X)  =  e“a"i(<0  1)LX 


(1 


The  final  complex  integral  is  evaluated  in  the  Appendix, 
where  the  desired  probability  is  derived  as  a  series  of  Marcum 
O-functions.  A  numerical  analysis  of  this  series  is  also  given 
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there,  on  the  basis  of  which  the  figures  of  the  present  section 

f 

have  been  produced.  When  there  is  no  signal  present,  the 
characteristic  function  simplifies  greatly,  and  the  contour 
(closed  in  the  lower  half  plane)  includes  only  the  simple  pole  at 
X  =  -i.  The  residue  at  this  pole  is  easily  evaluated,  with  the 
resulting  formula  for  the  probability  of  false  alarm 

PFA  =  (l/io)1'1  e~(*°~1)L  • 

For  a  single  pulse,  this  formula  and  the  corresponding 
formula  for  Pq  and  identical  to  the  expressions  for  PFA  and  PD 
for  conventional  detection  with  a  single  radar  hit.  This  is  to  be 
expected,  since  they  are  conditioned  on  given  loss  factors,  and 
the  CFAR  effect  associated  with  threshold  estimation  has  been 
eliminated  by  our  simplying  approximation. 

The  general  character  of  the  formulas  is  illustrated  in 
Figs.  2  through  7,  in  which  PD(a)  is  plotted  vs  a,  for  several 
sets  of  values  for  the  parameters  L  and  J.  In  all  cases,  the  PFA 
has  been  set  to  10“^.  The  dashed  curves  on  these  plots 
represent  the  PD  vs  SNR  for  a  conventional  detector  using  a 
single  radar  hit,  or  Marcum's  Q-function  for  N  =  1.  The  two 
curves  are  extremely  close  when  J  and  L  are  large,  and  when  the 
number  of  pulses  is  large,  they  are  relatively  insensitive  to  L. 
We  conclude  that  the  performance  of  the  ALR  test  can  be  expected 
to  be  very  similar  to  that  of  the  known-noise  test  of  Section  2, 
with  two  additional  losses:  a  CFAR  loss  typical  of  the  parameter 
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PROBABILITY  of  DETECTION  vs  SNR 


Probability  of  Detection  vs  SNR 


PROBABILITY  of  DETECTION  vs  SNR 


Probability  of  Detection  vs  SNR 


PROBABILITY  of  DETECTION  vs  SNR 


Probability  of  Detection  vs  SNR 


PROBABILITY  of  DETECTION  vs  SNR 


PROBABILITY  of  DETECTION  vs  SNR 
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L,  and  a  SNR  loss  which  is  essentially  that  predicted  by  the  Beta 
distribution  for  a  single  pulse. 
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(6)  SUMMARY 

It  has  been  the  intention  of  this  two-part  study  to 
discuss  the  problem  of  radar  operation  in  non-stationary 
interference,  from  the  point  of  view  of  final  target  detection. 
Those  aspects  of  the  radar  design  which  make  it  possible  to 
achieve  a  satisfactory  degree  of  interference  rejection  are  not 
analyzed  here,  but  are  taken  as  given.  We  refer  here  to  the 
provision  of  an  adequate  number  of  auxiliary  RF  channels  and  the 
means  to  produce  precisely  controlled  weighted  sums  of  their 
outputs.  The  choice  of  these  weights  is  often  discussed  in  terms 
of  interference  rejection,  or  nulling,  on  a  pulse  by  pulse 
basis.  By  viewing  the  problem  as  one  of  target  detection, 
utilizing  the  returns  from  a  sequence  of  pulses  which  form  a 
coherent  processing  interval  (CPI),  we  have  obtained  a  more 
complete  decision  algorithm.  This  algorithm  contains  both  the 
rule  for  choosing  the  auxiliary  weights  and  the  procedure  for 
combining  the  weighted  outputs  for  each  pulse  to  form  an 
integrated  resultant  for  target  declaration. 

Although  the  total  interference  is  modeled  as  Gaussian 
noise,  the  correlation  properties  of  this  interference  are 
presumed  to  be  unknown.  Detection  is  only  possible  if  these 
correlation  properties  can  be  estimated,  and  for  this  purpose  we 
have  assumed  the  availability  of  other  data,  namely  the  outputs 
of  adjacent  range  gates,  which  are  taken  to  be  signal  free.  The 
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important  additional  assumption  is  made  that  the  interference 
present  in  these  other  data  is  independent  of,  but  statistically 
identical  to  that  of  the  main  range  gate  in  which  target 
detection  is  attempted.  They  can  therefore  provide  a  valid  data 
base  from  which  the  covariance  matrices  of  the  interference  can 
be  estimated. 

In  this  analysis,  the  available  input  data  has  been 
characterized  statistically  and  the  methods  of  statistical 
decision  theory  have  been  applied  to  derive  a  detection 
procedure,  based  on  the  totality  of  original  inputs.  The 
resulting  procedure  can  be  considered  to  have  three  components, 
namely  covariance  estimation,  interference  rejection  and  coherent 
integration.  The  first  two  components  of  this  algorithm  are  in 
agreement  with  standard  procedures.  Covariance  estimation  is 
accomplished  by  means  of  the  sample  covariance  matrix  of  the 
secondary  (target-free)  data,  and  interference  rejection  is 
performed  by  sample  matrix  inversion  and  the  application  of  the 
corresponding  colored  noise  matched  filter  weights.  The  third 
component  of  the  decision  procedure  is  a  form  of  weighted 
coherent  integration  of  all  the  pulses  of  the  CPI,  in  which  the 
weights  are  dependent  on  estimates  of  the  residual  noise  level 
after  interference  rejection  on  each  individual  pulse.  It  was 
this  aspect  of  the  problem  that  was  of  most  interest  at  the 
outset,  although  it  is  reassuring  to  have  the  eminently 
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reasonable  and  conventional  noise  estimation  and  interference 
rejection  procedures  appear  as  derived  results  from  the  theory. 

In  Part  I  the  interference  was  allowed  to  have  arbitrary 
correlation,  not  only  among  the  multiple  RF  input  channels,  but 
also  among  all  the  pulses  of  the  CPI.  Because  of  the  great  number 
of  unknowns  this  entails,  the  number  of  secondary  data  inputs 
required  for  noise  estimation  must  be  very  great.  On  the  other 
hand,  due  to  the  generality  of  the  formulation,  an  exact  decision 
rule  was  obtained,  together  with  an  exact  evaluation  of  its 
performance.  The  latter  was  given  in  terms  of  probability  of 
detection  and  false  alarm,  as  functions  of  the  system 
parameters.  Since  no  special  structure  was  assumed  by  which  the 
data  from  separate  pulses  could  be  distinguished,  all  these  data 
were,  in  effect,  lumped  into  one  input  vector  in  which  a  target 
is  sought,  and  a  set  of  secondary  target-free  vectors  for  noise 
estimation.  Thus,  by  a  suitable  re-interpretation,  Part  I  can 
also  be  said  to  deal  with  the  problem  of  target  detection  based 
on  a  single  pulse. 

In  Part  II,  the  multiple  pulse  problem  is  specifically 
treated,  with  the  additional  and  reasonable  assumption  that  the 
interference  is  independent  from  pulse  to  pulse.  The  correlation 
properties  of  this  interference  are  otherwise  totally  unknown, 
and  must  be  estimated  from  the  target-free  inputs.  The  number  of 
these  secondary  inputs  needed  to  assure  adequate  performance  is 
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now  much  smaller,  being  related  to  the  number  of  RF  channels 

-f" 

themselves,  and  not  to  the  product  of  this  number  and  the  number 
of  pulses,  as  in, the  former  case.  On  the  other  hand,  the  analysis 
is  inherently  more  difficult,  and  two  distinct  approximations  had 
to  be  made  to  obtain  useful  results.  The  first  was  an 
approximation  to  the  decision  rule  itself,  which  seems  to  be 
justified  not  only  by  the  quantitative  arguments  made  in  its 
selection,  but  also  by  the  reasonable  character  of  the  resulting 
decision  rule,  the  approximate  likelihood  ratio  (ALR)  test.  The 
other  approximation  was  made  in  the  performance  analysis,  as  a 
result  of  which  these  results  are  somewhat  incomplete.  The 
limitations  of  this  analysis  were  discussed  in  section  5,  where 
the  insight  gained  from  Part  I  was  used  to  assess  their  probable 
impact . 

The  general  conclusions  of  the  study  are  best  expressed  in 
terms  of  the  performance  of  an  idealized  radar  which  uses  the 
same  primary  inputs,  but  has  the  advantage  of  knowing  the 
correlation  properties  of  the  total  noise.  This  radar  can 
therefore  dispense  with  the  secondary  inputs  of  the  real  system. 
The  decision  procedure  for  the  idealized  case,  derived  in  Section 
2,  consists  of  interference  rejection  on  each  pulse,  followed  by 
a  weighted  coherent  integration  of  the  pulses  of  the  CPI. 
Interference  rejection  is  based  on  the  known  noise  covariance 
matrices,  and  the  integrator  weights  are  inversely  proportional 
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to  the  residual  noise  levels  after  this  portion  of  the 
processing.  This  detector  could  be  described  as  a  conventional 
nulling  processor  applied  to  each  pulse,  followed  by  a  colored 
noise  matched  filter  which  performs  integration  over  the  CPI.  If 
the  external  interference  is  nulled  well  below  the  internal 
noise,  the  integrator  weights  will  be  equal,  and  the  second 
portion  of  the  processor  becomes  a  simple  integrator.  The 
advantage  of  using  weighted  integration  over  the  CPI,  compared  to 
the  use  of  constant  weights,  will  depend  on  the  residual  noise 
levels  and  their  variation  from  pulse  to  pulse.  No  general  rule 
can  be  given  for  the  performance  gain  in  this  case,  but  it  will 
generally  represent  an  improvement,  relative  to  the  use  of  an 
incorrect  set  of  weights. 

In  the  actual  radar  processor,  noise  estimates  are  used, 
both  in  the  nulling  of  individual  pulses  and  in  the  establishment 
of  the  final  detection  threshold.  The  use  of  estimated  covariance 
matrices  in  the  nulling  process  results  in  a  loss  of  signal  to 
noise  ratio  on  each  pulse.  This  loss  has  been  shown  to  be 
statistically  identical  to  that  derived  by  Reed,  Mallett  and 
Brennan,  in  a  well-known  paper<3).  The  relationship  between 
this  SNR  loss  and  the  number  of  secondary  inputs  (range  gates  in 
our  case)  is  well  understood,  on  the  basis  of  the  Beta 
distribution  to  which  these  losses  are  subject.  The  seco.d  use  of 
noise  estimation  relates  to  the  normalization  of  weights  for  the 
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coherent  integrator,  and  the  final  detection  threshold  itself. 

The  threshold  estimation  provides  an  approximately  CFAR  detection 
test,  and  also  leads  to  a  CFAR  loss  in  performance.  Because  of 
the  approximations  referred  to  already,  the  CFAR  loss  is  not 
accurately  assesed  by  our  results,  but  arguments  have  been  given 
that  it  should  not  be  large.  It  has  also  been  shown  that  this 
loss  should  be  approximated  by  the  loss  of  a  simple  linear  CFAR 
detector,  whose  threshold  is  based  on  K-M+l  noise  samples,  where 
K  is  the  number  of  secondary  inputs  and  M  the  number  of  RF 
channels  of  the  actual  system. 


APPENDIX 


EVALUATION  OF  THE  APPROXIMATE  DETECTION  PROBABILITY 


In  Section  5  an  expression  was  derived  for  the  approximate 
probability  of  detection,  conditioned  on  given  values  of  the  loss 
factors.  This  probability  was  expressed  as  the  contour  integral 


where 


oo— ie 

Prob(f>0)-PD(a)=gLj  *(X)^ 

— oo— ie 


$(X)  =  e-a 


i>LA  +  l5x 


The  detection  probability  depends  on  the  conditioning  only 
through  the  signal  parameter,  a,  defined  in  Section  5. 

For  a  >  0,  the  integrand  has  an  essential  singularity  when" 
1  equals  -i,  a  simple  pole  at  the  origin  and  a  pole  of  order  J-l 
at  i/(  2.o~l  )•  The  path  of  integration  is  completed  in  the  lower 
half  plane  and  then  shrunk  to  a  small  circle  about  the  essential 
singularity.  We  can  therefore  write 


1 

27Ti 


dA 

A 


where  the  contour  is  a  small  circle,  now  enclosing  the  point  -i 
in  a  positive  sense. 

We  make  the  definition 


y  53  (^o— i) L 

and  the  change  of  variable 

A  s— i(l— t)  . 
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which  brings  the  singularity  to  the  origin.  The  result  is 


„  /  ^  e~a~y  f  eyt  +  ?  dt 


|tl-e 

For  a  =  0,  the  singularity  at  the  origin  becomes  a  simple 
pole,  and  the  residue  easily  yields  the  formula  for  PFA  given  in 
Section  5.  In  general,  the  integral  will  not  yield  a  closed  form 
answer,  and  a  series  expansion  is  required.  The  difficulty  is 
caused  by  the  factor  exp(yt),  which  causes  another  essential 
singularity  at  infinity. 

We  make  the  following  expansion 


[<0-(<o-Dtr*  x*0/  m=0' 

whose  convergence  is  assured  by  the  fact  that  the  magnitude  of  t 
is  constant  and  arbitrarily  small  on  the  contour.  Then  we  find 


l-:I  |Vt+! 

2tti  J 

W-* 


The  contour  integrals  now  remaining  have  been  discussed  in 
detail  in  Ref.  1,  in  relation  to  the  Marcum  Q-function.  Let 

t 

N-l 


fN(u,a)  s  e  u  a  2  In-i(S^) 


and 
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FN(y,a)s  J  fN(u,a)du  _ 

7 

Then  Fn  is  equivalent  to  Marcum's  Q-function,  the  probability 
of  detection  after  non-coherent  integration  of  N  hits.  The 
parameter  y  is  the  normalized  threshold  and  a  is  the  normalized 
total  (i.e.  integrated)  signal  to  noise  ratio. 

In  Ref.  1  it  is  shown  that 

(  1  -  Fm+l(a-y)  :m^° 
l  F_m(y,a)  ;m<0 

where  the  evaluation  for  non-negative  m  has  been  made  by 
replacing  the  variable  of  integration,  t,  by  1/t.  Using  this 
result,  we  obtain 

P0(a)  =  ^y  ^(y.a) 

When  a  equals  zero,  this  formula  reduces  properly,  since 

FN(0,y)  =  l  , 

and 

FN(y,0)  =  e"7  . 

The  resulting  expression "for  the  probability  of  false  alarm. 


e~a~7  feyt  +  “  £dt  = 
2ni~  J  1-t 

|t|=«e 
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is  easliy  solved  for  £0,  given  PFA,  by  the  Newton-Raphson 
method . 

For  J  =  1,  the  detection  probability  is  simply 

PD(a)  =  Fi(y.a) 

since  all  the  binomial  coefficients  vanish  in  this  case  (no 
expansion  was  necessary  for  J  =  1).  This  formula  describes  the 
detection  performance  of  a  simple  coherent  integrator,  with  total 
SNR  equal  to  a  and  fixed  threshold,  and  we  see  that  the  CFAR  loss 
is.  not  accounted  for,  an  expected  consequence  of  the  use  of  our 
simplifying  approximation. 

For  other  values  of  J,  and  in  particular  for  large  J,  this 
formula  looks  quite  different,  since  will  be  a 

small  factor,  and  many  terms  of  the  series  will  contribute.  As 
the  results  presented  in  Section  5  show,  however,  our  formula 
produces  a  rather  conventional  PD  vs  SNR  curve  for  a  wide  range 
of  values  of  the  parameters,  L  and  J. 

Numerical  results  are  obtained  by  truncation  of  the  series 
and  the  use  of  a  truncation  bound,  which  will  now  be  derived.  We 
write  the  series  in  the  form 
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is  the  general  term  of  the  series  and  ffj  is  the  error  of 
truncation  at  the  term: 


We  now  make  use  of  the  identity 

F„(x.y)  =  FN_i(x>y)  +  fN(x,y)  , 

which  follows  easily  from  the  integral  representation  given 
above,  and  which  implies  that  the  sequence  FN(a,y)  increases 
monotonically  with  N,  since  the  PDF  ffl  is  necessarily 
non-negative.  The  factors 

[1  -  Fm(a,y)]  . 

then  decrease  monotonically,  and  we  shall  use  this  property  to 
write 


m 


After  replacing  the  index  m  by  M+m,  we  find  that  the 
inequality  becomes 


where  S  stands  for  the  series 


( J +M— 2)1  (M+m)!  \  l0  j 


But 


_  _  v  /j+M+m-2\  Mi  ml  /  \  m 

m  /  (M+m)l  \  l0  / 
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4-^) 
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—  ir 


and  hence 

eM  <  [1  -  Fu(a,y)](,+“~2)(<o-l)“  =  '0%  ■ 

The  truncation  bound  is  therefore  easily  computed  along  with  the 
general  term  of  the  series,  using  recursion  for  the  binomials  and 
simple  powers,  and  Shnidman's  algorithm ( 2)  for  the  Marcum 
functions.  This  procedure  has  been  used  to  produce  the  figures  of 
Section  5. 
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